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We study general properties of the Landau-gauge Gribov ghost form-factor o~(p 2 ) for SXJ(N C ) 
Yang-Mills theories in the d-dimensional case. We find a qualitatively different behavior for d = 3, 4 
with respect to the d = 2 case. In particular, considering any (sufficiently regular) gluon propagator 
T>{p 2 ) and the one-loop-corrected ghost propagator, we prove in the 2d case that the function c{p 2 ) 
blows up in the infrared limit p — > as — X>(0) ln(p 2 ). Thus, for d = 2, the no-pole condition 
a(p 2 ) < 1 (for p 2 > 0) can be satisfied only if the gluon propagator vanishes at zero momentum, 
that is, T>(0) = 0. On the contrary, in d = 3 and 4, u(p 2 ) is finite also if T>(0) > 0. The same results 
are obtained by evaluating the ghost propagator Q(p 2 ) explicitly at one loop, using fitting forms for 
T>(p 2 ) that describe well the numerical data of the gluon propagator in two, three and four space- 
time dimensions in the SU(2) case. These evaluations also show that, if one considers the coupling 
constant g 2 as a free parameter, the ghost propagator admits a one-parameter family of behaviors 
(labelled by g 2 ), in agreement with previous works by Boucaud et al. In this case the condition 
c(0) < 1 implies g 2 < g 2 , where g 2 is a "critical" value. Moreover, a free-like ghost propagator in 
the infrared limit is obtained for any value of g 2 smaller than g 2 , while for g 2 — g 2 one finds an 
infrared-enhanced ghost propagator. Finally, we analyze the Dyson-Schwinger equation for cr(p 2 ) 
and show that, for infrared-finite ghost-gluon vertices, one can bound the ghost form-factor a(p 2 )- 
Using these bounds we find again that only in the d = 2 case does one need to impose V(0) — 
in order to satisfy the no-pole condition. The d = 2 result is also supported by an analysis of the 
Dyson-Schwinger equation using a spectral representation for the ghost propagator. Thus, if the 
no-pole condition is imposed, solving the d — 2 Dyson-Schwinger equations cannot lead to a massive 
behavior for the gluon propagator. These results apply to any Gribov copy inside the so-called first 
Gribov horizon, i.e. the 2d result 23(0) = is not affected by Gribov noise. These findings are also 
in agreement with lattice data. 
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I. INTRODUCTION 

Green functions of Yang-Mills theories are gauge-dependent quantities. They can, however, be used as a starting 
point for the evaluation of hadronic observables (see for example [1-4]). Thus, the study of the infrared (IR) behavior of 
propagators and vertices is an important step in our understanding of QCD. In particular, the confinement mechanism 
for color charges [5] could reveal itself in the IR behavior of (some of) these Green functions. This IR behavior should 
also be relevant for the description of the deconfinement transition and of the deconfined phase of QCD. Indeed, at 
high temperature color charges are expected to be Debye-screened and the (electric and magnetic) screening masses 
should be related to the IR behavior of the gluon propagator (see for example [6-8]). 

Among the gauge-fixing conditions employed in studies of Yang-Mills Green functions, a very popular choice is 
the Landau gauge, which in momentum space reads p^A" (p) = 0. From the continuum perspective this gauge has 
various important properties, including its renormalizability, various associated nonrenormalization theorems [9] and 
a ghost-antighost symmetry [1]. In the past few years many analytic studies of Green functions in Landau gauge have 
focused on the solution of the Yang-Mills Dyson-Schwinger equations (DSEs), which are the exact quantum equations 
of motion of the theory (see for example [1, 2, 10, 11]). Since the DSEs are an infinite set of coupled equations, 
any attempt of solving them requires a truncation scheme. Then, some Green functions (usually the gluon and the 
ghost propagators) are obtained self-consistently from the considered equations, while all the other Green functions 
entering the equations are given as an input. For the coupled DSEs of gluon and ghost propagators two solutions have 
been extensively analyzed (see for example Chapter 10 in Ref. [5] and Rcf. [12] for recent short reviews). The scaling 
solution [1, 13-16] finds in d = 2, 3 and 4 an IR-enhanced ghost propagator G(p 2 ) and a vanishing gluon propagator 
V(p 2 ) at zero momentum. 1 In particular, the IR behavior of the two propagators should be given respectively by 
G(p 2 ) ~ (p 2 )-* ' 1 and by V(p 2 ) ~ ( p 2)2 KD +(2-d)/2 with Kq = Kd _ .2(d - 1) [14, 16]. On the other hand, the 
massive solution [18-26] gives (for d — 3 and 4) a free-like ghost propagator in the IR limit, i.e. hq = 0, and a massive 
behavior for the gluon propagator, 2 that is, T>(0) > and hd = (d — 2)/4. 

The existence of two different types of solution for the coupled gluon and ghost DSEs is now understood as due to 
the use of different auxiliary boundary conditions. 3 These conditions can be given in terms of the value of the ghost 
dressing function .F(p 2 ) = p 2 G(p 2 ) at a given momentum scale p [22, 31]. In particular, if one considers p = 0, it is 
clear that 1/J r (0) = gives an IR-enhanced ghost propagator G(p 2 ) while 1/.F(0) > yields a free-like behavior for 
G{p 2 ) at small momenta. As stressed in Ref. [22] (see also the discussion in Section 4.2.2 of Ref. [12]), the scaling 
condition 1/.F(0) = relies on a particular cancellation in the ghost DSE which, in turn, implies a specific "critical" 
value g 2 for the coupling constant g 2 [12, 32]. Thus, at least from the mathematical point of view, there is a one- 
parameter family of solutions for the gluon and ghost DSEs, labelled by g 2 or, equivalently, by 1/J"(0): in the case 
g 2 = g 2 one recovers the scaling solution while, for all cases g 2 < g 2 , the solution is a massive one. 4 In Ad 1 the SU(3) 
physical value of the coupling seems to select 5 the massive solution [21, 34]. 

At this point we should recall that, when considering gauge-dependent quantities in non-Abelian gauge theories, 
one has to deal with the existence of Gribov copies [35] (see also Ref. [36] for a recent review). Indeed, for compact 
non-Abelian Lie groups defined on the 4-sphere [37] or on the 4-torus [38] , it is impossible to find a continuous choice 
of one (and only one) connection A^{x) on each gauge orbit. The effect of Gribov copies is not seen in perturbation 
theory [35], i.e. the usual Faddeev- Popov-quantization procedure is correct at the perturbative level. However, these 
copies could be relevant at the non-perturbative level, i.e. in studies of the IR properties of Yang-Mills theories. 

Different approaches have been proposed in order to quantize a Yang-Mills theory while taking into account the 
existence of Gribov copies (see for example [35, 39-45]). The one usually considered, both in the continuum and on the 



1 For the explanation of color confinement based on the scaling solution see for example Section 3.4 in [17] and references therein. 

2 The possible existence of a dynamical mass for the gluons, as well as its relation to quark confinement through vortex condensation, has 
been discussed a long time ago in Ref. [27] . 

3 The possibility of having different non-perturbative solutions for DSEs in relation with different boundary conditions has been discussed 
for example in Refs. [28, 29] (see also Section 1.4 in [30] and Sections 3 and 3.1 in Ref. [1]). 

4 If g 2 > g 2 one gets a negative ghost propagator at small momenta. 

5 Recently, in Ref. [33] , it has also been shown — using a renormalization-group approach — that in three and four space-time dimensions 
only the decoupling solution is expected to be physically realized. 
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lattice, is based on restricting the functional integration to a subspace of the hyperplane of transverse configurations. 
The original proposal, made by Gribov [35], was based on the observation that the Landau gauge condition d fl A b {x) = 
allows for (infinitesimally) gauge-equivalent configurations if the Landau-gauge Faddeev-Popov operator M bc (x, y) — 
—5(x — y)d^D b ^ has zero modes. (Here D b £ is the usual covariant derivative.) Indeed, since an infinitesimal gauge 
transformation Sco(x) gives A^x) — > A^x) + D b foj c (x), it is clear that the exclusion (in the path-integral measure) 
of the zero modes of A4 (x, y) implies that gauge copies connected by such infinitesimal gauge transformations are 
ignored in the computation of expectation values. In order to exclude these zero modes, Gribov considered a stronger 
condition by requiring that the functional integration be restricted to the region O of gauge configurations A^x) 
defined as 

fi = { Al{x) : d^(x) = , M bc (x, y) > 0} . (1) 

This set, known as the (first) Gribov region, clearly includes the vacuum configuration A^x) = 0, for which the 
Faddeev-Popov operator is given by —S(x — y)S bc d 2 . The region Q can also be defined (see for example [46, 47]) as 
the whole set of local minima 6 of the functional £ [A] = J d d x A b ^ (x) A b ^ (x) . Since usually each orbit allows for more 
than one local minimum of £[A], it is clear that the region fi is not free of Gribov copies. On the contrary, in the 
interior of the so-called fundamental modular region A, given by the set of the absolute minima of the functional £[A], 
no Gribov copies occur [43, 48]. 

The characterization of the fundamental modular region A, i.e. finding the absolute minima of the energy functional 
£[A], is a problem similar to the determination of the ground state of a spin glass system [49, 50]. Thus, a local 
formulation of a Yang-Mills theory, with the functional measure delimited to A, is not available, whereas a practical 
way of restricting the physical configuration space to the region Q was introduced by Gribov [35]. To this end, he 
required that the ghost dressing function F{p 2 ) cannot have a pole at finite nonzero momenta. After setting 

p z p z 1 — a(jr) 

this condition can be written as 

a(p 2 ) < 1 for p 2 > , (3) 
where <r(p 2 ) is the so-called Gribov ghost form-factor [35]. Indeed, since the ghost propagator is given by 

p) , (4) 

i.e. it is related to the inverse of the Faddeev-Popov matrix M bc (x, y), the above inequality — known as the no-pole 
condition — should be equivalent to the restriction of the functional integration to the Gribov region f2, i.e. to the 
condition M bc (x,y) > 0. 

From the discussion above, it is clear that both scaling and massive solutions of DSEs satisfy the no-pole condition, 
i.e. 1/F{p 2 ) = 1 — cr(p 2 ) > for p 2 > 0. Indeed, in the scaling case [1, 13-16], this condition [together with the 
condition 1/.F(0) = 0] is imposed from the beginning to the solution of the DSEs. On the contrary, for the massive 
solution, the no-pole condition is either verified a posteriori, as in Ref. [20], or used [together with the condition 
1/J"(0) > 0] as an input for the solution of the DSEs, as in Refs. [19, 25, 26]. In particular, in Ref. [26], the value of 
1/F(0) is fixed using lattice data. 

The restriction to the first Gribov region f2 is also always implemented in lattice numerical simulations of Green 
functions in Landau gauge by (numerically) finding local minima of the functional £ [A] . Results obtained using very 
large lattice volumes [51-54] (see also Chapter 10 in Ref. [5], Section 3 in Ref. [12] and Ref. [55] for recent short 



G(P 2 ) 



N 2 - 1 



, be 



For this reason this gauge condition is often indicated in the literature as minimal Landau gauge. 
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reviews) have shown that in d — 3 and 4 the gluon propagator T>(p 2 ) is finite and nonzero in the limit p — > while 
the ghost propagator G{p 2 ) behaves as 1/p 2 - On the contrary, for d = 2 the lattice data [56-59] are in quantitative 
agreement with the scaling solution and one finds n D = k g 0.2. 

Since the region f2 is not free of Gribov copies, their (possible) influence on the numerical evaluation of gluon 
and ghost propagators has been studied by various groups [61-66]. It has been found that these effects are usually 
observable only the IR limit and that any attempt to restrict the functional integration to the fundamental modular 
region A gives a stronger suppression at small momenta for both propagators, i.e. reducing the value of T>(0) and 
increasing that of 1/J r (0). More recently, it has been suggested [67-70] that the one-parameter family of solutions 
obtained for the gluon and ghost DSEs should be related 7 to Gribov-copy effects and that the value of 1/J r (0) could 
be used as a gauge-fixing parameter. This analysis finds indeed IR-cnhanccd ghost propagators (and sometimes 
a disconcerting over-scaling 8 ). On the other hand, the gluon propagator still shows a finite nonzero value at zero 
momentum, that is, T>(0) > 0. Moreover, this approach does not explain why the numerical results found in d = 2 arc 
different from those obtained in d — 3 and 4, even though Gribov copies inside the first Gribov region f2 are clearly 
present in any space-time dimension d > 1. 

From the analytical point of view, following Gribov's approach, Zwanziger modified the usual Yang-Mills action 
in order to restrict the path integral to the first Gribov region ft [71] Although this restriction is obtained using a 
non-local term, the Gribov-Zwanziger (GZ) action 9 can be written as a local action and it is proven [72-74] to be 
renormalizable. At tree level the GZ gluon propagator is given by T>(p 2 ) — p 2 /(p 4 + A 4 ), where A is a parameter 
with mass-dimension 1. At the same time, the ghost propagator is given by G(p 2 ) ~ 1/p 4 - Thus, as in the scaling 
solution of the gluon and ghost DSEs, the gluon propagator is null at zero momentum 10 and the ghost propagator 
is IR-enhanced [42]. These tree-level results, also in agreement with the original work by Gribov [35], have been 
confirmed by one- loop calculations in the three- and four-dimensional cases [77-81]. 

More recently, the GZ action has been modified by considering (for d — 3 and 4) dimension-two condensates [82- 
85]. The corresponding action, called the Refined Gribov-Zwanziger (RGZ) action, still imposes the restriction of the 
functional integration to the region and it is renormalizable. However, the RGZ action allows for a finite nonzero 
value of T>(0) and for a free-like ghost propagator G(p 2 ) in the IR limit. Thus, nonzero values for these dimension-two 
condensates yield for the gluon and ghost propagators an IR behavior in agreement with the massive solution of the 
gluon and ghost DSEs. 11 Indeed, the RGZ tree-level gluon propagator describes well the numerical data in the SU(2) 
case [59, 60], for d — 3 and 4, and in the SU(3) case [88] with d — 4. It is also interesting to note that the fitting values 
for the dimension-two condensates are very similar for the SU(2) and SU(3) gauge groups in the four-dimensional 
case. 

As stressed above, the restriction of the functional integration to the first Gribov region and the no-pole condition 
(3) are key ingredients in the study of the IR sector of Yang-Mills theories in Landau gauge. However, to our knowledge, 
a detailed investigation of the properties of the Gribov form- factor a(p 2 ) as well as of the possible implications of the 
no-pole condition was missing up to now, although some interesting one-loop results were already presented in Refs. 
[84, 89-91]. In particular, in Appendix B.2 of [84] it was shown that, if the gluon propagator T>(p 2 ) is positive, then 
in the 2d case the derivative da(p 2 )/dp 2 is negative for all values of p 2 , i.e. <r{p 2 ) is largest at p 2 — 0. Also, in Ref. 
[89] it was proven that in the RGZ framework the form-factor a(p 2 ) presents a logarithmic IR singularity — ln(p 2 ) for 
d = 2. This result precluded the use of the RGZ action in the two-dimensional case, leading to a first interpretation 
of the different behavior found in lattice numerical simulations for the 2d case, compared to the d = 3 and 4 cases. 
Similar findings have been (more recently) presented in Refs. [33, 92, 93]. 

In this work we collect some general properties of the Gribov form- factor (r(p 2 ) and we study the consequences 
of imposing the no-pole condition. In particular, in Section II, using the expression for o{p 2 ) obtained from the 



7 This identification is, however, based on several (unproven) hypotheses, as already stressed in Ref. [58]. 

8 Let us recall that the scaling solution is supposed to be unique [31] . 

9 See again Ref. [36] for a comprehensive review of the GZ action. 

10 One can, however, obtain a finite nonzero value for T>(0) within the GZ approach by considering a non-analytic behavior for the free 
energy of the system [75, 76]. 

11 Let us note here that a massive behavior for these two propagators has also been obtained in Refs. [86, 87] using different analytic 
approaches. 
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FIG. 1. Feynman diagrams for the one-loop-corrected Landau-gauge ghost propagator. Dashed lines represent ghosts, the curly 
line represents gluons. 

evaluation of the ghost propagator at one loop, we prove that <j(p 2 ) attains its maximum value at p 2 = for any 
dimension d > 2. Since this expression for a(p 2 ) depends on the gluon propagator T>(p 2 ), in the same section we also 
investigate (for a general <i-dimensional space-time) which IR behavior of the gluon propagator is necessary in order to 
satisfy the no-pole condition cr(p 2 ) < 1. By considering a generic (and sufficiently regular) gluon propagator D(p 2 ), we 
find in the d = 2 case that cr(p 2 ) is unbounded unless the gluon propagator is null at zero momentum. More exactly, 
we find a(p 2 ) ~ — 2?(0)ln(p 2 ) in the p — » limit, in agreement with [89]. This result does not apply to the d = 3 
and 4 cases. Indeed, in these cases one can introduce, for all values of p 2 , simple finite upper bounds for the Gribov 
form-factor. In Section III we present explicit one-loop calculations for cr(p 2 ) using for the gluon propagator T>(p 2 ) 
linear combinations of Yukawa-like propagators (with real and/or complex-conjugate poles), which have been recently 
used to model lattice data of the gluon propagator in the SU(2) case [59, 60]. Besides confirming the results obtained 
in Section II, we also find that the ghost propagator admits a one-parameter family of behaviors [21] labelled by the 
coupling constant g 2 , considered as a free parameter. Moreover, the massive solution Q(p ) ~ 1/p 2 , corresponding 
to cr(0) < 1, is obtained for all values of g 2 smaller than a "critical" value g 2 . At the "critical" value g 2 , implying 
<r(0) = 1, one finds an IR-enhanced ghost propagator. (As already stressed above, the case g 2 > g 2 corresponds to 
<t(0) > 1 and one obtains a negative ghost propagator at small momenta.) Finally, in Section IV, we analyze the 
DSE for <r(p 2 ). We stress that in this case we do not try to solve the DSE but we focus only on general properties of 
this equation. As we will see, considering IR-finite ghost-gluon vertices, we confirm and extend the one-loop analysis 
of the no-pole condition presented in Section II. In particular, after introducing bounds for the Gribov form-factor, 
we show again for d = 2 that the gluon propagator T>(p 2 ) must vanish at zero momentum in order to keep o(p 2 ) 
finite. On the contrary, such a constraint does not apply in the three- and four-dimensional cases. We also present 
alternative evidence for the d = 2 result using a spectral representation for the ghost propagator in the DSE. 

It is important to note that all our results in Sections II and IV apply irrespective of which set of Gribov copies 
(inside the region f2) is considered, i.e. they are not affected by the so-called Gribov noise. We end with our Conclusion 
in Section V. Some technical details have been collected in four Appendices. In particular, in Appendix B we present 
properties of the Gauss hypergeometric function 2-F1 (a, b; c; z) that are relevant to prove some of our results. 

II. THE ONE-LOOP-CORRECTED GHOST PROPAGATOR AND THE GRIBOV FORM-FACTOR 

In this Section, as well as in Section III below, we consider the one-loop-corrected Landau-gauge ghost propagator, 
diagrammatically represented in Figure 1. This propagator can be written [for the SU(A C ) gauge group in the 
e?-dimensional case] as 

where 8 ab T>{q 2 )V^ v {q) is the tree- level gluon propagator [not necessarily given by T>(q) = l/q 2 ] and V^q) = 
\Piiu — Q^Qv/l 2 ) is the usual projector onto the transverse sub-space, i.e. q^'P tlu (q) — 0. We have also considered the 
tree- level ghost-gluon vertex igf adc p V: where p is the outgoing ghost momentum. The indices a, d, c refer, respectively, 
to the incoming ghost, to the gluon and to the outcoming ghost. After using f adc f cdb = —N c 5 ab , valid for the adjoint 
representation, we obtain 



g( P 2 ) = - [i + <j(p 2 )} 



(6) 
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where a(p 2 ) is the momentum-dependent function 



a(p 2 ) = g 2 N< 



VnVv f d d g 1 
p 2 J (2ny(p-q) 2 



v{ q 2 )v^(q) ■ 



Finally, we can write [as in Eq. (2)] 



G(P 2 ) = 



1 



p 2 1 — u(p 2 ) 



(7) 



(8) 



which corresponds to the usual resummation of an infinite set of diagrams into the self-energy. Note that this 
resummation makes sense only when cr(fc 2 ) < 1, i.e. when the no-pole condition (3) is satisfied. 

Clearly the function a(p 2 ) is dimensionless and it should go to zero for p — > oo, modulo possible logarithmic 
corrections. Also, this function coincides with the so-called Gribov ghost form-factor [35, 90, 91], even though the 
latter is obtained in a slightly different way. 12 As discussed in the Introduction, the no-pole condition <j{p 2 ) < 1 for 
p 2 > should be equivalent to the restriction of the path integral to the first Gribov region f2 [defined in Eq. (1)]. In 
this section we will derive general properties of cr(p 2 ) in d > 2 space-time dimensions. In particular, as we will see 
below, the weaker condition a(p 2 ) < +oo is already sufficient to obtain a strong constraint on the IR behavior of the 
gluon propagator V{p 2 ) in the d = 2 case. 

A. First Derivative of a(p 2 ) for d — 2 



Following Appendix B.2 of Ref. [84] one can show that, if the gluon propagator T>(q 2 ) is positive, then in the 2d 
case and for all values of p 2 we have 



dp 2 



< 



(9) 



with <j(p 2 ) defined in Eq. (7). This implies that a(p 2 ) attains its maximum value at p 2 = 0. To this end, we choose 
the positive x direction parallel to the external momentum p and write (using polar coordinates) 



g 2 N c p 2 J (2n) 2 (p-q) 



7 V(q 2 )V^(q) 



qdq 
4tt 2 



V{q 2 ) 



2tt 



d0 



1 - cos 2 (<9) 



p2 _|_ q2 



2 qp cos(#) 



(10) 



The integral in d9 can be evaluated using contour integration on the unit circle and the residue theorem. This yields 



/ d9 

Jo 



1 - cos 2 {9) 



p 2 + q 2 — 2qp cos(#) 



2 - z 2 - z 2 



i dz 



4pq z 2 - z(p/q + q/p) + 1 p 2 



9(p 2 -q 2 ) + ^6(q 2 -p 2 ), (11) 



where § dz represents the integral on the unit circle \z\ = 1, we indicated with z the complex-conjugate of z = e l6 and 
9{x) is the step function. This integral is also evaluated in Eqs. (A17) and (A18) in the Appendix A (for the general 
d-dimensional case). Considering also Eq. (B8) and Eq. (A5), with d = 2, wc find 



<P 2 ) 
g 2 N c 



1 

47T 
1 

8^ 



-/ 



Jp 



dq 



o P 



V(q 2 ) 



P dr? 2 2 , f°° dq 2 n . 2 , 
o P 

dq 2 V(q 2 ) 



j q 

9{p 2 -q 2 ) , 9{q 2 -p 2 ) 



(12) 
(13) 
(14) 



12 Note that in the Gribov ghost form- factor there is usually an extra factor l/(Af? — 1) [90, 91] compared to our Eq. (7). However, this is due 
to the fact that in Eq. (5) above we considered for the Landau-gauge gluon propagator the usual expression D^Kq 2 ) = & ah V l _i V (q)'D{q 2 ) 
while, in the derivation of the Gribov ghost form-factor, one usually writes q) = w(g as in Eq. (255) of Ref. [90]. 
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Then, by using d x 9(x) = S(x), where S(x) is the Dirac delta function, the derivative of a{p 2 ) with respect to p 2 yields 

*W _ r d?2 ^r) _ _ g% f Af v(q2) _ (15) 



dp 2 8tt J ^ ' p A 8tt P 4 

which is clearly negative, for any value of p 2 , if T>(q 2 ) is positive. We can evaluate the limit p 2 — > of this derivative 
using, for example, the trapezoidal rule 13 

lim ^ = - Hm |^ t- [V(p 2 ) + 2,(0)] = - hm . (17) 

One arrives at the same result after writing Eq. (15) as 

W£± = - ^ t dx lW) = - ^ ^ b{i)) (18) 
9p 2 8V J ax U{xv > Sup 2 p 2 ' K > 

where D(p 2 ) is a primitive of T>(p 2 ), that is, D'(p 2 ) = V(p 2 ) where we indicate with ' the first derivative with respect 
to the variable p 2 . The limit p 2 — > then yields again Eq. (17). 

Clearly, one finds an IR singularity at p 2 — 0, unless T>(0) = 0. If this condition is satisfied, using again the 
trapezoidal rule, we have from Eq. (17) that 

lim = - lim ^ £- V(p 2 ) = - lim ^ ^ = « = - ^ lim V W . (19) 

p^o dp 2 P *->o 8tt P 4 2 v ; P ^o 16tt p 2 16tt p *-x> y ' v ; 

For a gluon propagator V(p 2 ) that is regular at small momenta, i.e. that can be expanded as V{p 2 ) ss V(0)p 2 + 
T>" '(0)p 4 /2 at small p 2 , the above limit is finite. On the other hand, if the leading IR behavior of T>(p 2 ) is proportional 
to p 2v with 1 > rj > 0, as found for example in the 2d case in Refs. [14, 16, 56, 59], then the above limit gives a 
singular value, due to the non-integer-power (and non-analytic) behavior of V(p 2 ). 



B. Infrared Singularity of a(p 2 ) for d = 2 



Here we prove that — for d — 2 and for any gluon propagator T)(p 2 ) that goes to zero sufficiently fast at large 
momenta, e.g. as 1/p 2 , and that is reasonably regular at small momenta, e.g. that can be expanded at p — as 
V{p 2 ) w D(0) + Bp 2 ^ + Cp 2!i (with £ > rj > and V{0),B and C finite) 14 — the ghost form-factor (7) displays a 
logarithmic divergence for p — »■ proportional to ©(0). Indeed, by considering Eq. (12), one obtains 



9 2 N C 



1 


■ /-p 


47T 


Vo 


1 


lim 


87T A^co 


1 


lim 


on A— >oo 


1 


lim 


Sir A^oo 



S & { jf 7 + 2 jf 7 c(0) + 2 jf 7 1 B< ' 2) - B,0, l 

f p2 da; 

/ — V(x) + V(0) In 
Jo P 

£(p 2 ) - £(0) 



T v{q2) 



A 2 



+ 2 



dg 

V g 2? ? + M 



V{q 2 ) - V(0) 



2-q 



(q 2r > + M) 



+ ^(0) In ( ^ 



dx ■ 



,.7,-1 

x r > + M 



V{x) - D(0) 



{x 71 + M) 



(20) 
(21) 
(22) 
(23) 



The trapezoidal rule gives the numerical approximation 

f dx f(x) = [ f(b) + f(a) ] + 0(b - a) 3 , (16) 

which can be obtained by integrating f(x) rj /(a) + (x — a) [/(b) — f(a)] /{b — a). Thus, the trapezoidal rule is equivalent to using a 
linear Taylor expansion f(a) + (x — a)f'(a) for f(x) with the first derivative /'(a) approximated by a (first forward) finite difference 
[m -f(a)]/(b-a). 

For example, in Eq. (116) below, the Taylor expansion of T>(p 2 ) at p 2 = is of the type considered here with § = 1. 
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where x — q 2 , D(x) is again a primitive of T>(x) and M > is a (finite) constant. If we indicate with H{x) the 
quantity in square brackets in the last line, then we have 



/ ^ 2, -^°i + P ( ) ln(^-'^)ln(^ + M ) A " 



- / dx ln(x" + M) ff'(x) 

n Jv 2 



(24) 



Note that for V(x) = l/(x T > + M) we find H(x) = 1/M = V(0) and the last term in Eq. (24) is zero. Since 
Hindoo V(x) = we also have that lim^oo H(x) = T>(0) and the two logarithmic singularities for infinite A cancel 
each other. Thus, we get 



a(p 2 



l 2 N r , 8ir 



-I D{lp2) 2 D{0) - V(0) In (p 2 ) + I H{p 2 ) In (p 2 " + M) 



P 



n 



+ 



i r 

V Jp2 



dx In {x r > + M) 



rjM [D(x) - X>(0)] - x(xi + M) V{x) 
3.1+77 



(25) 



If V{x) <~ 1/a; at large x, it is easy to check 15 that a(p 2 ) is null for p 2 — > 00, as expected. At the same time, in the 
limit p 2 — > we obtain 



a(0) 



1 



f 



5 2 A/ C 8tt 



— <^ £>(0) - lim D(0) ln(p 2 ) + - i? (0) ln(M) 



da; In {x 71 + M) 



r]M[V(x) - X>(0)] - x (x" + M) ©'(a;) 



,1+7) 



(26) 



where we used 



lim £^l_m = &{ 0) = V(0) 
p 2 -)-0 p 2 



(27) 



and H(0) = —MB is a finite constant. 16 Finally, in Appendix C we show that, under the assumptions made for 
the gluon propagator, 17 the last term on the r.h.s. of Eq. (26) is finite. Thus, the only IR singularity in the ghost 
form-factor a{p 2 ) is proportional to —T>(0) ln(p 2 ). This result is in qualitative agreement with [89]. An IR singularity 
plaguing the 2d calculation has also been recently obtained in Ref. [92]. 

An alternative (equivalent) proof 18 can be done by performing an integration by parts. Then, Eq. (12) becomes 



<r(p 2 ) _ 1 
g 2 N c 8tt 

1 

~~ 8^ 
1 

~ 8tt 



p dx f 00 dx 

J p2 X 



P 

Dip 2 ) ~ D(0) 
p 2 

D(p 2 ) - D(0) 



+ In (a;) V(x) 



dx ln(x) V'{x) 



- \n(p 2 ) V{p 2 ) - dx In (a) V'{x) 
Jp 2 



(28) 
(29) 
(30) 



15 See details in Appendix C. 

16 Note that here we used the IR expansion T>(p 2 ) ss T>(0) + Bp 2v + Cp 2 ^ for the gluon propagator. 

17 In the same Appendix we will also show that the hypotheses considered above for the gluon propagator T>(x) can be relaxed. 

18 Note that both proofs are singling out the singularity —T>(0) ln(p 2 ) by essentially doing a Taylor expansion of T>(p 2 ) at p 2 = 0. 



9 



where we used the assumption T)(x) <~ 1/x at large x. Note that the above result coincides with Eq. (25) when M = 0, 
which implies = 2?(0) — V(x). A second integration by parts yields 



9 2 N C 



ln(p 2 ) V{p 2 ) - [x ln(x) - x] V'(x) 



dx [x In (a;) - x] V"(x) 



(31) 



d(p 2 ) - m 

J_\ D(p 2 ) - D(0) 
8n\ p 2 

Here we used the hypothesis that T>'(x) goes to zero sufficiently fast at large momenta, e.g. as 1/x 2 . As before, one 
easily sees that <r(p 2 ) is null for p 2 — »■ oo (see Appendix C). At the same time, under the assumptions made for the 
gluon propagator T>(p 2 ), in the limit p 2 — > we obtain 



\n(p 2 ) V(p 2 ) + [p 2 \n(p 2 ) - p 2 ] V'{p 2 ) + dx[x ln(x) - x] V"{x) \ . (32) 



g 2 N c 



87T ^ p 2 ^0 J 



(33) 



and we again find 9 an IR singularity proportional to — X>(0) ln(p 2 ), unless one has V(0) = 0. 



Thus, in the 2d case and using a generic (sufficiently regular) gluon propagator, a null value for 2?(0) is a necessary 
condition to obtain a finite value for cr(0) at one loop. As a consequence, the condition T>(0) = must be imposed if 
one wants to satisfy the no-pole condition (3) and keep the functional integration inside the first Gribov region il. It 
is important to stress again that our proofs apply to any Gribov copy inside the first Gribov horizon, i.e. the result 
T>(0) — is not affected by the Gribov noise. 



C. Properties of a(p 2 ) in d Dimensions: Approximate Calculation 



We can easily extend the result 



da( P 2 ) 
dp 2 



< 



(34) 



to the d-dimensional case by using for the integral in d d q the so-called y-max approximation or angular approximation 
(see for example [13, 18, 94]). The same approach allows us to show that the IR singularity —T>(0) ln(p 2 ) is present 
only in the two-dimensional case. Indeed, by using hyperspherical coordinates (see Appendix A) and by considering 
the positive X\ direction parallel to the external momentum p, we can write the c?-dimensional ghost form-factor (7) 
as 

<?(P 2 ) _ V»Vv f d d q 1 ,^ _ r^q 4 ' 1 ^^ [ jc. 1 -™* 2 ^) 



(35) 



g 2 N c p 2 J {2n) d (p-q) 2 ^ >' ^ J 1 (2n) d yH ' J (p - q) 2 ' 

In the y-max approximation one substitutes l/(p — q) 2 with 1/p 2 , for q 2 < p 2 , and with 1/g 2 , for p 2 < q 2 . Then, we 
obtain 



<y(p 2 ) i 

g 2 N c {2v) d 



dq 



■V(q 2 ) + 



p 2 V1 ' ' J P ~* q 2 



dq- 



■V(q 2 ) 



/[ 



l- cos 2 (</»!)] dn d 



n d d-i 



(27r) d 2d 

n d d-i 

(27r) d 2d 



dq 



2 q 



■V{q 2 



d-2 

dq 2 ^V(q 2 ) 

q 



dq 2 q d - 2 V(q z ) 



0(P 2 ~ q 2 ) Q(q 2 -P 2 ) 



(36) 
(37) 
(38) 



In Appendix C we will prove that the integral on the r.h.s. of Eq. (33) is finite under the hypotheses made for the gluon propagator 
V(p 2 ). In the same Appendix we will also show how these hypotheses can be relaxed in this case. 
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where we used Eq. (A12). Note that, for d — 2 and using Eq. (A5) one gets the exact result (14). By repeating the 
argument shown in the Section II A, the proof of the inequality (34) follows directly from Eq. (38). 
At the same time, we can write Eq. (37) as 



a(p 2 ) Q d d-l 



J* dqq^^l + J°° dqq d - 3 V(q 2 ) 



Up 2 ,™), (39) 



g 2 N c [2w) d d 

where the integral I d (p 2 ,£) is defined in Eq. (B34). In Appendix B we have also shown that, for d > 2, this integral 
is finite when the gluon propagator V(p 2 ) is finite and nonzero at p 2 = 0. Thus, using the y-max approximation, 
we find that only in the 2d case the condition T>(0) = is necessary in order to obtain a finite value for the Gribov 
form- factor a(p 2 ) for all values of p 2 . 

Of course, in case of ultraviolet (UV) divergences we should regularize the integral defining <r(p 2 ), as done for 
example in Section III C below for the Ad case using the modified minimal subtraction (MS) scheme and dimensional 
regularization. One can also consider a fixed momentum \x and subtract 20 the value cr(/i 2 ) from the Gribov form- factor 
<j{p 2 ). Due to the use of the y-max approximation the result of the subtraction is very simple. Indeed, instead of Eq. 
(39) we have the relation 



a(p 2 ) - a{n 2 ) n d d-l 



fP 2 x d/2-l ft* 2 rn 2 x d l 2 ~ X 

\ dx — V(x) + / dxx d/2 - 2 V(x) - / da; V(x) 

Jo P 2 Jpi Jo M 



g 2 N c ~ [2n) d 2d 
which is valid for p 2 < fi 2 as well as for [i 2 < p 2 . Then, we find again 

da(p 2 ) _ 2at Q d d-l f , x^ 2 - 1 



(40) 



1 fP r d/2-L 

- / dx t—V(x) < (41) 

Jo P 



dp 2 9 c (2n) d 2d 

if V(x) is positive. We can also easily check that, for V(0) > and d > 2, the Gribov form-factor a(p 2 ) in Eq. (40) 
does not display an IR singularity. 

D. Properties of a(p 2 ) in d Dimensions: Exact Calculation 

One can improve the results obtained in the previous Section by considering the formulae reported in Appendices 
A and B, which allow us to perform the angular integration in Eq. (35) without approximations. Indeed, we have 21 

^=r**WT* V ^ / dild 2 + V C °: 2(0l) rn - ^ M. oo) , (42) 

9 2 N C Jo (2TT) d J p A + q A - 2pq cos(0i) 

with I(p 2 ,v,d,£) defined in Eq. (B31). Since v = 1 in this case, for d > 2 we can also make use of the inequalities 
(B36) and write 

I d (p 2 ,™) < /(p 2 ,l,d,oo) < I d (p 2 ,™) . (43) 



2 (d-l) 

Note that I d (p 2 , oo) is the same integral obtained on the right-hand side of Eq. (39). Thus, the y-max approximation 
of the previous Section provides, for d = 3 and 4, an upper bound for the Gribov ghost factor. On the contrary, 
for d = 2, the above inequalities become equalities. At the same time, as one can see in Appendix B, the integral 
I d (p 2 , oo) is finite (for d > 2) also if T>(0) is nonzero, i.e. we do not need to impose the condition T>(0) — in order 
to attain a finite value for a(p 2 ) in the IR limit. 



20 This is equivalent to a momentum-subtraction (MOM) rcnormalization scheme defined by the condition G(fJ- 2 ) = l/^ 2 - 

21 Again we make the hypothesis that a regularization is introduced in the case of UV divergences. In Appendix D we explicitly show how 
to extend the proof to a MOM scheme, i.e. by subtracting the value ct(/x 2 ) from a(p 2 ) where fj, is a fixed momentum. 
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By evaluating the derivative with respect to p 2 of the result (B32) we also obtain 



1 da (p 2 ) 



g 2 N c Op 2 



(2tt) c 



0(p 2 - q 2 ) 



2 F 1 (l,l-d/2;l + d/2;q 2 /p 2 ) 



+ 



q 2 0(p 2 -q 2 ) 
p 6 

8(q 2 -p 2 ) 



2 F{(l,l-d/2;l + d/2;q 2 /p 2 ) 



2 F[(l,l-d/2;l + d/2-p 2 /q 2 ) 



(44) 



where 2 F[ (a, 6; c; z) indicates the derivative with respect to the variable z of the Gauss hypergeometric function 
2-Fi (a, b; c; z) (see Appendix B). Here we used again the properties of the theta and of the Dirac delta functions and 
Eq. (B6). For d = 2, the last two terms in Eq. (44) are null [see Eq. (B23)] and, using the result (B8), we find again Eq. 
(15). In the Ad case one can use the expression (B9) for the Gauss hypergeometric function 2 F\ (1, 1 — d/2; 1 + d/2; z). 
Then, from Eq. (44) — or, equivalently, by evaluating the derivative with respect to p 2 of Eq. (Dl) in Appendix D 
— we find that 



da(p 2 ) _ g 2 N c 

Op 2 32 7T 2 

9 2 N C 
32 it 2 



dqV(q 2 ) 



2„3 



2<f - 3p 2 q 

p 6 



dq 



dy V{y 2 p 2 ) (2y 5 - 3y 3 



P(<? 2 ) 

q 

V{q 2 



dq 



(45) 



(46) 



where y = q/p and we have used Eq. (A5). For T>(p 2 ) > both terms in square brackets are negative, i.e. the 
derivative d<j{p 2 ) / dp 2 is negative for all values of the momentum p. Let us note that in the original work by Gribov 
[35] the same result was proven [see comment after Eq. (37) in the same reference] under the much stronger hypothesis 
of a gluon propagator V{q 2 ) decreasing monotonically with q 2 over the main range of integration. 

A similar analysis can be done in the 3d case using Eq. (B17). In order to simplify the notation we define 



*(*) = 2*1 (1,-1/2; 5/2; z) 



3 

8z 



V5 



arcsinh f 



VV 1 - z 



This gives 



16z 3 



z (z — 3) + \fz (3 — 2z — z 2 ) arcsinh 



(47) 



(48) 



Then, after setting d — 3 in Eq. (44) and using Eq. (A5), we obtain 



da(p 2 ) = g 2 N c 
dp 2 



3tt 2 



dqq 2 V(q 2 



0(P 2 ~ q 2 ) 
P 4 



q 



g 2 N c 

6ir 2 



r 2 



+ 



P z 



p< 

0(q 2 - P 2 ) 



q 2 8(p 2 -q 2 ) 
p 6 



p° 



f°° da; n . , T . 



(49) 



(50) 



where we also made the substitution x = q 2 . Next, the change of variable x — yp 2 in the first integral and x = p 2 /y 
in the second integral yield 

d<j(p 2 ) _ g 2 N c j 



dp 2 



6tt 2 



dy^V(yp 2 ) 
P 



f ^^-V(p 2 /y)V{y)\. 
Jo P Vv J 



(51) 



As one can see in Figure 2, the factor — [ ^(y) + y W(y) ] is negative for y e [0, 1] . At the same time, from Eq. (B24) 
we know that *&'(y) is negative for y > (see also the corresponding plot in Figure 2). Thus, for a positive gluon 
propagator T>(p 2 ), the 3d derivative da(p 2 )/dp 2 is negative for p 2 > 0. 
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0.2 0.4 0.6 0.8 



FIG. 2. The functions -[*(?/) + yV'{y)] (full line) and (dashed line) for y £ [0, 1]. 



Finally, we can consider a general d > 2 and, after suitable changes of variables (for p 2 > 0), we write 



1 da(p 2 ) n d d-lp 



d-i ( A 



dy y 2 - d ' 2 V{p 2 /y) 2 F{ (1,1- d/2; 1 + d/2; y) 



5 2 iV c dp 2 (27r) d d 2 

- dytf" 2 - 1 TUvp*) 2 Fj. ( ] . "J -d/2: I - ,/■•.>:,,! F, /2 F;(l.l-d/2:l + ( //2://) J> . (52) 



Note that the dependence on p 2 is only in the global factor p d ~ 4 and in the argument of the gluon propagator. From 
Appendix B we know that the derivative iF[ (1, 1 — d/2; 1 + d/2; x) is negative for x G [0, 1] and d > 2 and that, 
under the same hypotheses, the expression in square brackets is positive. Thus, for a positive gluon propagator, both 
term in the r.h.s. of the above expression are negative and we have proven that, for any dimension d > 2, the Gribov 
form-factor cr(p 2 ) (at one loop) is monotonically decreasing with p 2 , i.e. it gets its maximum value at p 2 = 0. 



III. EVALUATION OF THE ONE-LOOP CORRECTED GHOST PROPAGATOR USING (LINEAR 
COMBINATIONS OF) YUKAWA-LIKE GLUON PROPAGATORS 

In the previous Section we have proven that, at one- loop level and for a sufficiently regular gluon propagator D(p 2 ) 1 
the Gribov ghost form- factor (j(p 2 ) is always finite in three and four space-time dimensions while, in d = 2, one needs 
to impose T>(0) = in order to avoid an IR singularity of the type — T>(0) ln(p 2 ). In this Section we present an explicit 
calculation of a(p 2 ) at one loop for d = 2, 3 and 4 using, for the gluon propagator, results recently presented in Ref. 
[59, 60] from fits to lattice data of T>(p 2 ) in the SU(2) case. The expressions obtained below for the ghost propagator 
G(p 2 ) will be used in a subsequent work [95] to model lattice data of SU(2) ghost propagators. 

In this Section, besides recovering the same results reported in Section II, we also find that the ghost propagator 
G(p 2 ) admits a one-parameter family of behaviors [21, 34] labelled by the coupling constant g 2 , considered as a free 
parameter. The no-pole condition cr(0) < 1 implies g 2 < g 2 , where g 2 is a "critical" value. Moreover, for g 2 smaller 
than g 2 one has <t(0) < 1 and the ghost propagator shows a free-like behavior in the IR limit, in agreement with the 
so-called massive solution of gluon and ghost DSEs [18-25]. On the contrary, for g 2 = g 2 one finds cr(0) = 1 and the 
ghost propagator is IR enhanced [1, 13-16]. 
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A. Yukawa-Like Gluon Propagators and Set Up 



In Ref. [59, 60] the SU(2) gluon propagator was fitted in 2, 3 and 4 space-time dimensions using, respectively, the 
functions 22 



V{p 2 ) = C 



p 2 + IjP + s 

p 4 + u 2 p 2 +t 2 ' 



(53) 



V(p 2 ) = C 



p 4 + (s + l)p 2 + s 



p 6 + (k + u 2 )p A + (ku 2 + t 2 )p 2 + kt 2 



and 



V{p 2 ) = C 



p 2 + s 



p 4 + u 2 p 2 + t 2 



(54) 



(55) 



The last two propagators are tree-level gluon propagators that arise in the study of the RGZ action [82-85]. The 
first one is a simple generalization of the form (55) that fits well the 2d data. Note that these three functions can be 
written as a linear combination of propagators of the type 1/ (p 2 + oj 2 ), where lu 2 is in general a complex number. [In 
the 2d case we need to consider the more general form p v /(p 2 + w 2 ) with r\ > 0.] Thus, in order to evaluate a(p 2 ) in 
Eq. (7) using the above gluon propagators T>(p 2 ), we first consider the integral 



V^Vv f d d g 1 

p 2 J (2n) d (p - q) 2 q 2 



.2 \ 



Q 2 



(56) 



The evaluation of f(p,ui 2 ) can be done in three and four space-time dimensions by introducing Feynman parameters 
and applying the usual shift in the momentum q. The integration then yields 



1 



{^) d > 2 Jo 
1 



dx 



A d/2 - 2 T(2-d/2) 



f dx [ x dy [le d/2 - 2 r(2 - d/2) + x 2 P 2 e d/2 - 3 r(3 - d/2) 

Jo Jo 1.2 



(57) 



with 



A = -x 2 p 2 + xp 2 + (1 - x)uj 2 



©22, 2 i 2 

= —x p + xp +yoj 



(58) 



Since the Gamma function has the behavior T(x) rj l/i for small x, it is clear that the first two integrals are UV 
finite for d < 4 while the third one is UV finite for d < 6. Below we will calculate the integral (57) for d = 3 and 4. 
We start from the case d = 3, where all terms are finite, and then we evaluate the integral for the case d = 4, using 
the MS scheme. On the contrary, as stressed above, in the 2d case one needs to evaluate the more general function 



f(P,u 2 ,v) 



P^lPv 



d 2 q 



1 



p 2 J (27r) d (p - q) 2 q 2 + u 2 



(59) 



with t) > 0. This case will be treated (in a slightly different way) in Section HID. 

Most of the analytic results reported in this Section have been checked using Mathematica and/or Maple. 



Note that, for consistency with the notation used in Ref. [59], in this Section the non-integer power of the momentum p is r\ and not 2r] 
as in the rest of the manuscript. 
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B. Ghost Propagator in the 3d Case 



In the 3d case the residual x- and y-intcgrations in Eq. (57) are straightforward and give 



r 1 ( 


v y 




- — arctan 


+ 


47rp \ 







(p 2 - ur)Vw 



+ 



7rp 



{p 2 + u 2 ) 2 j p 

arctan 



32np 2 u) 2 64ttw 2 32np 3 u) 2 



3(p 2 -w 2 )^ 3p A -2p 2 uj 2 + 3(u 2 ) 2 ( p \ 3wp 

H — — — n arctan 



32np 2 uj 2 



32irp 3 uj 2 



Uttlo 2 



(60) 



where the three square brackets highlight the contribution from the three terms in Eq. (57). Here we have only made 
the assumption p 2 > 0. By simplifying the above result, we find 



f(p 2 ,u 2 ) = — !— arctan 



4np 



P 



(p 2 — uj 2 )\/ui 
16irp 2 LU 2 



Tip 



2\2 



(p 2 ~u 2 ) 

32ttuj 2 16ttp 3 oj 2 



■ arctan 



P 



1 



32irp 3 uj 2 



g(p 2 ,^) 



where 



g{p\io 2 ) = -WP 4 + 2p 3 V^ - 2p(u; 2 ) 3 / 2 + 2 (p 2 + lo 2 ) 2 arctan ( -^=) 

In order to use the result (61) we need to write the gluon propagator (54) as 

a (3 7 



V(p 2 ) = 



P 2 +Uj\ p 2 + U)\ p 2 + UJ 3 



2 ' 



(61) 



(62) 



(63) 



Here u 2 , u 2 and wf are the three roots of the cubic equation, with respect to the variable p 2 , obtained by setting 
equal to zero the denominator of Eq. (54). Thus, by combining Eqs. (7), (63) and (56) we can write for the function 
<j{p 2 ) in the 3d case the expression 



a(p 2 ) = g 2 N c [af(p 2 ,uj 2 ) + Pf(p 2 ,co 2 ) + jf(p 2 ,u 2 " 



(64) 



or 



a(p 2 ) = 



9 2 N C 



32 7T wV 



[ag(p 2 ,uij) + (3g(p 2 ,uj 2 ) + 7 .9(p 2 ,w|)] 



(65) 



with g(p , uj 2 ) given in Eq. (62). In general, the roots uj 2 , lo 2 and uj 2 are all real or there is one real root, for example 
uj 2 , and two complex-conjugate roots, i.e. (w|)* = wf, implying also ft = 7*. Since the fits in Ref. [59, 60] support the 
latter case we write 



(3 = a + ib , 7 = a — ib 



(66) 



and 



uj 2 . = v + iw 



Wo = v — iw . 



Then, following for example [96], we find for ui\ the relations 



2 _ 



y/v + iw = —= V \J v 2 + w 2 + v H — \/ \J v 2 + w 2 — v 

V2 V \/2 V 

(w|) 3/2 = (v + iw) 3/2 = (v + iw)Vv + iw , 
p p p 



/ujI Vv + iw \Jv 2 + w 2 



\V — ii 



(67) 



(68) 

(69) 
(70) 



15 



and similar results for wf . We also use the expression (see for example [97]) 



arctan (z) = — arg 



i + z 



.in 



I — z 



i + z 



Vz ^ {i, -i} 



This allows us to write the function a(p 2 ) only in term of real quantities, i.e. 



a(p 2 ) 



4 7T 0j\ p 3 



where g(p 2 ,u> 2 ) is given in Eq. (62) above. Also, we have 



f R ( P 2 ) = fi( P 2 ) + f2( P 2 ) + h{ P 2 ) + Up 2 ) + h( P 2 



with 



MP 2 ) 

f2{ P 2 ) 

h(p 2 ) 



av + bw 
~ P 2R 2 ' 

(av + bw) V-R + v — (bv — aw) y/R—v 



P 



1 ay/R + v 

~2 



by/R^ 



U(p 2 ) = A(p 2 ) 



V27T 

p* (av + bw) + 2 ap 2 R 2 + R 2 (c 



bw) 



2irR 2 p 3 



h(p 2 ) = -L(p 2 ) 



2 p A (bv - aw) + 2 bp 2 R 2 + R 2 (bv + aw) 



2nR 2 p 3 



and 



A(p 2 ) 



L(p 2 ) = In 



arctan (^p) 

7T + arctan ( ^[^ ) 

yV + 2p 2 v + R 2 



if R - p 2 > 
if R - p 2 < 



R + p (p + V2~VR-v) 



(71) 

(72) 

(73) 

(74) 
(75) 
(76) 
(77) 
(78) 

(79) 

(80) 
(81) 



R = \Jv 2 + w 2 . 

One can check that a(p 2 ) is null in the limit p — » oo. Finally, by expanding <r(p 2 ) around p 2 = in Eqs. (72)-(81) 
we find 

a(p 2 ) a °" d2 



g 2 N c 6tta 



^—9aR 2 - (av -bw) (2v-R) 



+VR-v 



24:V2wR 3 

9bR 2 - (bv + aw) (2v + R) a R 2 + 2 (av + bw)uj 2 



P + O (p 2 ) , 



(82) 



24V2^i? 3 Z2u 2 R 2 

which implies G(p 2 ) oc p~ 2 at very small momenta. However, if the constant term in the above expression is equal to 
l/(g 2 N c ), yielding <r(0) = 1, then one gets in the IR limit Q(p 2 ) oc p~ 4 or G(p 2 ) cx p~ 3 , depending on whether the 
term aR 2 + 2 (av + bw)oj 2 vanishes or not. In particular, in the original GZ case, i.e. when the terms containing 
ui\ and a are absent, we do recover the usual 1/p 4 behavior. Also note that for purely imaginary poles, i.e. when 
v = b = (and R = w), the condition <r(0) = 1 simplifies to 



V2a 



67T 



(83) 
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Clearly, for a given value of N c and with a suitable choice of g 2 , one can always set c(0) = 1 in Eq. (82). For 
example, using the numerical data in the second column of Table XI of Ref. [59] and iV c = 2, we find from Eq. (82) 
the result 23 

^Q- w 0.039(0.001) - 0.017(0.003) p + O (p 2 ) (84) 

and we have <r(0) = 1 if g 2 w 12.82. Thus, if one considers g 2 as a free parameter, then Eq. (82) gives a one-parameter 
family of behaviors, labeled by g 2 . For a specific value of g 2 — g 2 we have a(0) = 1 and one finds an IR-enhanced 
ghost propagator at one loop. On the contrary, for g 2 < g 2 we obtain <r(0) < 1 and G(p 2 ) oc p~ 2 in the IR limit. 
Finally, for g 2 > g 2 the no-pole condition <r(0) < 1 is not satisfied, i.e. the ghost propagator is negative in the IR 
limit. These findings are in qualitative agreement with the DSE results obtained in Refs. [21, 34]. Finally, note that 
at small momenta the function a(p 2 ) in the above formula (84) is decreasing as p 2 increases, as expected from Section 
II D. 



C. Ghost Propagator in the 4d Case 



We want now to evaluate f(p 2 ,uj 2 ) in Eq. (57) for d = 4. As stressed above, in this case we have to deal with UV 
divergences. We do the calculation in the MS renormalization scheme using dimensional regularization with d = 4 — e. 
For the first term in Eq. (57) we have 



(47T) 



2-d/2 



167T 2 



dx 



A d / 2 - 2 r(2-d/2) 



16» 2 Jo 



dx 



- 7_e + In (4tt) - In (A) 



where je is the Euler constant. Then, using the usual MS prescription, we find 



-1 



Jo 



dx In 



x 2 p 2 + xp 2 + (1 - x)uj 2 



16tt 2 

and the da; integration yields 



-1 
16tt 2 



mli 



16^ 2 J 



dx 





--)] 


(x-\ 


P 2 )_ 



(85) 



(86) 



1 + In 1 + 



-1 

16p 2 7T 2 



-2p 2 + p 2 In 



+ oj 2 In 



p 2 +L0 2 



(87) 



where \x is the renormalization scale. For the second term in Eq. (57), which is also divergent, we first perform the y 
integration exactly, obtaining 



r(2-d/2) 



(47T) d / 2 CJ 2 (d-2) Jo 

The £ expansion then gives 
1 



/ dx ( y —x 2 p 2 + xp 2 ) 
Jo 



2 \d/2-l 



32tt 2 



dx (1 — x) 



In 



-x 2 p 2 + xp 2 
~~P 2 



xp 



xp z 



In 1 



d/2-l 



1 



xp" 



where we have already applied the MS prescription, and after integrating in dx we find 



1 



192jAj 2 tt 2 



p 4 (p 2 + 3c 2 ) In ( ^ ) + {p 2 +u 2 f\n 



(p 2 + LU 2 \ 


2 2 


{ U 2 ) 


+ p UJ 







-7p'-uj z + 3p 2 \n[^ 



(88) 



(89) 



(90) 



The errors in brackets have been evaluated using a Monte Carlo analysis with 10000 samples (see Ref. [59] for details). 
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Finally, the third term, which is finite, yields 



111 ( X + ) - 111 (*) 



1 



96p 4 w 2 7r 2 



yr'^V - + 2(/+c, 6 )ln(^±^ 



(91) 
(92) 



By summing the three results above we ultimately find (in the MS scheme) 

1 



64p 4 W 2 7T 2 



[h(p 2 ,uj) + f 2 (p 2 ,u 2 ) + / 3 (pV 2 )] 



with 



f 1 (p 2 ,u;)=p 4 (u; 2 -p 2 )ln(^j , 

Mp 2 , lo) = - (/ - P W + 3 p 2 cu 4 + c 6 ) In 

h{p 2 ,u) = p 2 u 2 



p 2 +UJ 2 



5 P 2 + lu 2 +p 2 ln( ^ ) -4p 2 In 



p 2 +UJ 2 



(93) 

(94) 
(95) 
(96) 



As shown in Ref. [59, 60], in the Ad case the fit of the gluon-propagator data is done using the expression (55). Thus, 
in order to use the above result (93)-(96), we need to write the gluon propagator as 



V{p 2 ) 



p 2 + W 2 p 2 + LU 2 _ 



(97) 



where co± are the roots of the quadratic equation, with respect to the variable p 2 , obtained by setting equal to zero 
the denominator of Eq. (55). Then, the ghost form- factor in the MS scheme is given by 



a™(p 2 ) = g 2 N c [ a+ f{p 2 ,u\) + a_f{p 2 ^ 2 _)~ 



and we have 



1 - <J MS {p 2 ) 



(98) 



(99) 



Note that the function a ms (p 2 ) is real. From [59, 60] we know that bj\ are complex-conjugate roots, i.e. w 2 = {oj\) 
and a- = a* + . By writing a± — a ± ib and uj± = v ±iw we find 



<^{p 2 ) = [-P 2 h( P 2 ) + R%(p 2 ) + P -%(p 2 ) p -%(p 2 )] 



with 



h(p 2 ) = (av + bw)[£ 2 {p 2 ) + £ 3 (p 2 )} - (bv - aw)[ ai {p 2 ) - a 2 {p 2 )] , 
t 2 (p 2 )=a[5 + £ 1 (p 2 ) + £ 2 (p 2 )+£ 3 (p 2 )-4£ i (p 2 )} - b[ ai (p 2 ) - a 2 (p 2 ) - 4a 3 (p 2 )} , 
hip 2 ) = [1 - 3£ 3 (p 2 )}(av 3 - bwv 2 + vaw 2 - bw 3 ) - 3a 2 (p 2 )(bv 3 + awv 2 + vbw 2 + aw 3 ) , 
t 4 (p 2 ) = £ 3 (p 2 )(av 4 - 2wbv 3 - 2vbw 3 -aw 4 ) + a 2 (p 2 )(bv 4 + 2awv 3 + 2vaw 3 - bw 4 ) 



(100) 

(101) 
(102) 
(103) 
(104) 
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and 

„2 



flip 2 ) = In (JjJ , (105) 
< 2 (P 3 ) = In (J;) , (106) 
+ (107) 

^)=h,( ^ + y + g ). (los) 

ai (p 2 ) = arctan , (109) 

wp 2 
i? 2 + i>p 2 1 
w 

^v+p 2 



n 2 (p 2 ) = arctan f — ) . CLIO) 
a 3 (p 2 ) = arctan ( - W 2 ) , (111) 



R = Vv 2 + w 2 . (112) 
Also note that, at large momenta, one gets 



•■v'-^K?)- (ll3) 

Finally, by expanding <r MS (p 2 ) around p 2 = in Eqs. (100)-(112) we obtain 

a ~Ms(p2^ 6 a In (=l) — 6 6 arctan (^) — 5 a 
5 2 iV c ~~ 64^ 

-ll + 61n^^ (av + wb) + 6(bv - aw) arctan (f) 



+ ■ 



IM^JP (I") 
Thus, if cr MS (0) = 1 we have that G MS (p 2 ) ~ 1/p 4 at small momenta (plus logarithmic corrections). Clearly, also in 
4d, we obtain a one- parameter family of behaviors, labelled by the value of g 2 , and the IR-enhanced ghost propagator 
corresponds to the upper value of g 2 allowed by the no-pole condition (3). With the numerical values reported in the 
second column of Table IV of Rcf. [59] and N c = 2 we find 24 

a MS( 2\ 

K1 2 ' = 0.0240(0.0007) + [-0.0082(0.0003) + 0.0060(0.0002) In (p 2 )] p 2 (115) 
2 9 

and the condition cr MS (0) = 1 corresponds 25 to g 2 s=s 20.83. Note again the negative sign of the leading order corrections 
at small momenta (see Section II D). 

D. Ghost Propagator in the 2d Case 

As stressed in Section III A above, Ref. [59] has shown that the fit of the gluon-propagator data in the 2d case can 
be done using the expression 

v{p } = ~p r Tu\ + ~p^TJT ' (116) 



24 Again, the error in parentheses have been evaluated using a Monte Carlo analysis with 10000 samples. 

25 Clearly, different renormalization schemes will modify the constant term in Eq. (115) and the value of g 2 . 
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where c is real, a_ = a*,, uj 2 _ = (uji_)* and u>± are the roots of the quadratic equation, with respect to the variable 
p 2 , obtained by setting equal to zero the denominator of Eq. (53). Thus, in order to evaluate the ghost form-factor 
<j(p 2 ) we need to consider the function f(p,uj 2 ,ri), defined in Eq. (59) above. To this end, we can choose again the 
positive x direction parallel to the external momentum p and consider polar coordinates. Then, after evaluating the 
angular integral we find 



f(p,u 2 ,v) = 



47T 



f p dq g 1+r > r°° dq _L 

Jo P 2 q 2 +u 2 J q 1 - r iq 2 +to 2 



(117) 



valid both for 77 = and for 77 > 0. 

In the case 77 = the momentum integration is straightforward giving 



f(p,0J 2 ) = lim 



1 



A->oo 4ir 

1 



f 

Jo 



dq q 



p 2 q 2 + uj 2 



1 



dq 



q q 2 + lu 2 



= lim 



A-s-oo Air \ 2p 2 



In 1 + 



1 

8^ 



1 

P 2 



In 1 + 



+ 



+ 



In 1 + 



In (A) - In (p) - 1 In (A 2 + u 2 ) + 1 In (p 2 + u 2 ) 



p 2 



(118) 

(119) 
(120) 



Note that the second term above blows up logarithmically in the IR limit p — > 0, in agreement with the result obtained 
in Section II B. 

For X) > the second integral in Eq. (117) can be written, after the change of variable t = uj 2 / (q 2 + uj 2 ) , &s 



dq 



1 



B 



OJ 



1 



where 



' p q 1 -^q 2 +UJ 2 2{lo 2 Y-^/ 2 \p 2 +UJ 2, 2' 2 



B(x;a,b) = f dtf 1 - 1 (l-t)^ 1 
Jo 



(121) 



(122) 



is the incomplete Beta function, which is defined for a, b > [97], implying 2 > 77 > in our case. For the first 
integral in Eq. (117) we cannot use directly the changes of variable v = l/q and t = 1/(1 + uj 2 v 2 ) because we get an 
incomplete Beta function (122) with b < 0. In this case it is convenient to introduce a Feynman parameter (using 
non-integer exponents) and write 



f p dq _q*_ 1 = J_ / _ rn / 

J p 2 q 2 -^ q 2 +uJ 2 p 2 \ 2J J 

2 



dxx-"' 2 



q 3 dq 



[q 2 + (1 - x)uJ 2 f- V ' 2 



(123) 



[p 2 + (l-x)oj 2 ] 1 - r ' /2 



+ - [p 2 + (1 - x)co 
rj 



2 11/2 - [(1 - x)uj 2 ] v/2 



(124) 



where we have also done the integration in dq. After suitable changes of variables, the last formula can be written as 



dq q 3 



1 



1 



p 2 q 2 V q 2 + uj 2 



B 



77 77 

2( w 2)l-rj/2 " \ p 2 + UJ 2' 2 ' 2 



+ 



p 2 + lu 2 



rip 2 (w 2 ) 1 "^ 



B 



UJ 



-1-5 i + H) - 

P 2 +uj 2 ' 2' + 2 ] 



r>W)^ B l}-l l + \) " (125) 
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where B (a, b) = B (1; a, b) = T(a) T(b)/T(a + b) is the Beta function. Thus, by summing the two results above, we 
find 



(a; 2 )"/ 2 l p 2 +m 2 

47T?7P 2 UJ 2 
( w 2)„/2-l 



47T 



r] V 2' 2/ 



2 i ,,,2' 1 O' T 9 



)-B( 1 -f, 1 + f) 



(126) 



+ 



(cj2) W2 p 2 + ^2 



4irrjp 2 lu 2 



(127) 



Note that for p = the incomplete Beta function B (w 2 /(p 2 + w 2 ); 1 —77/2,1 + 77/2) becomes the Beta function 
B (l — |, 1 + I). Then, by Taylor expanding f(p 7 uj 2 , 77) for small momenta p, we obtain 



yielding a constant contribution at p = 0. 

Using the expression (116), the ghost form-factor in the 2d case is given by 

a(p 2 ) = g 2 N c [a + f(p 2 ,L0 2 + ) + a.f{p 2 ,uj 2 _) + icf(p 2 ,L0 2 + , V ) - icf(p 2 ,L0 2 _, V )] , 



(128) 



(129) 



with f(p,u>) and f(p ,u> , 77) defined, respectively, in Eqs. (120) and (127). Of course, the function a{p ) is real. By 
writing a± = a±ib and to 2 . — v ±iw we get for the first two terms above 



1 I 1 



a + f(p 2 ,u 2 + ) + a.f{p 2 ,u; 2 _) = — [a£ 3 (p 2 ) + &a 2 (p 2 )] 



+ ^2 [( av + bw) £ 5 (p 2 ) - (bv - aw) a 3 (p 2 )] 



where l^{p 2 ), a 2 (p 2 ), a 3 (p 2 ) and R have already been defined in Eqs. (107), (110), (111) and (112) and 



h(p 2 ) = In 



p 2 



(130) 



(131) 



As shown in Section II B, there is a logarithmic singularity l^(p 2 ) ~ — ln(p 2 ) at small momenta proportional to the 
gluon propagator at zero momentum, that is, 2?(0) = 2(av + bw)/R 2 . We also have 



icf{p 2 ,u)\ ,rj) - icf(p 2 ,uj 2 _,T]) = -2cQ? [f{p 2 ,u\ ,77) 



(132) 



where we have indicated with 9 the imaginary part of the expression in square brackets. 

One can easily check that cr(p 2 ) is null at large momenta. Finally, the results (130) and (132), together with the 
expressions (128) and (129), allow us to evaluate the behavior of the ghost propagator at small momenta. We obtain 



a(p 2 ) 1 f ap 



bw 



g 2 N c 8tt [ 2i? 2 
-2cQ? 



i In I 4 



bv 



R 2 



R 2 



arctan 



47T77 



B 



1 + i 



p 1 



1 J ap 2 av + bw 
8^[2R^ + ' 



I hi ( 4 



—2c sin 



R 2 

(I - l) arctan (H) 



4^77(1 + 77/2)^ 
few — aw 



+ O (p 2+I ' 



+ o(p 4 ) 



(133) 



i? 2 



arctan 



/ w \ wp 



2cwp v 



47T 77(1 + 77/2) i? 2 



0(p 2+I ') . (134) 
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Note that, if <r(0) = 1, one finds a ghost propagator with a behavior l/p 2+v in the IR limit. As in 3d and in 4d we 
have a one parameter family of solutions labelled by the value of g 2 . 

As explained in Ref. [59], the 2d data for the gluon propagator suggest the relations a — —b and v = w, implying 
av + bw = and R 2 = 2v 2 . Then, we find 



°M - a* - 2 -° Id - Jl ^SF B - ?■ 1 + 3) - ^ife - s& + <p2+,) ■ ,135) 

Using the approximate result 7/ « 1 (see again Ref. [59] ) this formula simplifies to 



(2 v 2)„/2-l ^ ^ cp , ap 2 



_ ^ + 4^1-1/^2 ^ 
5 27V C _ 32« 6™ 32^ + U U> i ' ^ 

On the contrary, for N c — 2 and with the numerical values reported in [59] — see the second column of Table XIV 
and, for the exponent 77, the last line of Table XIII — we find for Eq. (134) the numerical results 

^Q- w 0.029(0.004) - 0.029(0.005) p°- 909(0 049) - 0.023(0.004) p 2 . (137) 
2.9 

The coefficient (av + bw)/R 2 cx 2?(0), multiplying the logarithmic IR singularity, is zero within error and we have 
omitted the corresponding term. Note that a(p 2 ) decreases for increasing momenta p 2 , as proven in Section II A 
above. Also note that we have <r(0) = 1 for g 2 w 17.24 and in this case the ghost propagator behaves as <~ 1/p 2 - 9 in 
the IR limit. 



IV. THE GHOST PROPAGATOR BEYOND PERTURBATION THEORY 

The one-loop analysis above has shown that, in the 2d case, an IR singularity —T>(0) ln(p 2 ) appears in the Gribov 
form-factor a(p 2 ) when p 2 — > 0. Thus, one needs a null gluon propagator at zero momentum in order to satisfy the 
no-pole condition <r(0) < 1. On the contrary, for d = 3 and 4, we found that a(p 2 ) is finite also for V(p 2 ) > 0. 

In this section we improve our analysis by considering the DSE for the ghost propagator G(p 2 ) (see for example 
[13, 18, 94]). As stressed in the Introduction, here we do not try to solve the ghost propagator DSE, but instead 
we concentrate on general properties of this equation for different space-time dimensions. In particular, the results 
obtained in Section II are confirmed by considering a generic (sufficiently regular) gluon propagator V(p 2 ) and an 
IR-finitc ghost-gluon vertex igf adc p\T\ u (p 1 q). 



A. The 2d Case 

In the 2d Landau gauge the DSE for the ghost propagator is written as 

^ = p 2 g 2 N c J ^ P xT )w (p, q )s fl V( q 2 )P^( q )g(s 2 ) , (138) 

where s = p — q, the gluon and the ghost propagators — respectively V(p 2 ) and Q(p 2 ) — are full propagators and we 
indicated with ig f adc p\T \ v (p, q) the full ghost-gluon vertex. The above result implies 

*{P 2 ) = ^ / ^PxTx v Ms v V{ q 2 )P, v { q )^ Y -^- ) (139) 

if one uses Eq. (8). For a tree- level ghost-gluon vertex T\ v (p, q ) = 5\ v and using the transversality of the gluon 
propagator we finally find 

<V) - «>^ p -f J ("») 
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which should be compared to the one- loop result (7). As in Section II A above, we can choose the x direction along 
the external momentum p obtaining (using polar coordinates) 

g 2 N c J Air 2 J Q s 2 [1 - a{s 2 )] 

with s 2 = p 2 + q 2 -2pq cos(0). 

This equation will be analyzed below using two different approaches. A first result can, however, be easily obtained 
using again the y-max approximation, as in Section II C above. This gives us 

g 2 N c Stt | J p 2 [1 - a(p 2 )} J p2 x[l- o{x)] J V > 

where we have done the angular integration and set x = q 2 . In the limit of small momenta p 2 we then obtain 

_ i r lim £ m±m + r dl . ri ^> . i , (143) 



g 2 N c 8tt {p*-x> 2 p 2 [1 - a(p 2 )] J x [1 - a(x)] 

In order to avoid IR singularities in the above equation we have to impose T>(p 2 ) ss Bp 2v , i.e. the gluon propagator 
should be null at zero momentum. In particular, if er(0) < 1, i.e. for a free-like ghost propagator at small momenta, 
it is sufficient to have r\ > 0. On the contrary, if the ghost propagator is IR enhanced and 1 — cr(0) cx x K for small x 
with k > 0, then the condition rj > n should be satisfied. Note that the predictions of the scaling solution [14-16], 
i.e. rj = 0.4 and k — 0.2, are consistent with the above inequality. The same results can also be obtained by setting 
p 2 = directly in Eq. (141). This makes the 9 integral trivial and gives 

(7(0) rqdq V{q 2 ) 



-s 

JO 



g 2 N c J 47T q 2 [1 - a(q 2 )} ' (W4) 

Note, however, that in both cases we essentially miss the logarithmic IR singularity — ln(p 2 ) which is found below. In 
the first case this is probably related to the very crude y-m&x approximation. On the contrary, in Eq. (144), this is 
due to the (improper) exchange of the q integration with the p 2 — > limit [89] . 



1. Bounds on the Gribov Form-Factor 



Since the Gribov form-factor is non-negative, we can easily construct a lower bound for the l.h.s. of Eq. (141) by 
writing 



°(P 2 ) 
g 2 N c 



> 



qdq 
4tt 2 



V(q 2 



d0 



1- cos 2 (<9) 



J 2 (p 2 ,oo) 



1 

4ir 



Jo P 2 



7(p 2 ,l,2,cx) 
dq 



■V(q 2 ) 



(145) 



%^ )+ r 

v Jp 1 

where we use the definitions (B31), (B34) and the relations (B36). The last integral in the above equation has already 
been analyzed in Section II B, where it was shown that ■, oo) develops an IR singularity proportional to — ln(p 2 ) 
if V(0) ^ 0. Thus, a(p 2 ) also is IR singular, unless X>(0) = 0. 

One can also find an upper bound for a(p 2 ) and check that the IR singularity is indeed only logarithmic. To this 
end we can notice that, if er(0) < 1, one can write 26 



<P 2 ) 
g 2 N c 



< 



Qdq 
4tt 2 



V{q 2 ) f 
Jo 



de- 



i 



s 2 W 
CT (o); 



^2(P 2 ,00) 

l-a(0) ' 



(146) 



Recall that, in the 2d case and in the one-loop approximation, the function c(p 2 ) is decreasing as p 2 increases, i.e. the maximum value 
of a(p ) is obtained for p = (sec Section II A). However, the proof presented here can be easily modified for the case when a(p ) < 1 
for all momenta p and the maximum value of c(p 2 ) is not attained at p = 0. Finally, one should note that in the DSE (140) one uses 
explicitly Eq. (8). Thus, when estimating the integral in Eq. (141), we cannot simply impose u(p 2 ) < +oo but we have to consider the 
stronger condition cr{p 2 ) < 1. 
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where we have also used Eqs. (145) and (145) above. Therefore, the upper bound also blows up as — ln(p 2 ) in the IR 
limit. At the same time, if cr(0) = 1, with <r(p 2 ) «1- cp 2K at small momenta we find 



°(P 2 ) r°°Qdq„, 2 , A 2 " 1 - cos 2 (0) 



f OO J fZTV 



„2+2k 



1-(t(s 2 ) 



(147) 



Note that the quantity in square brackets in the last integral is finite at s = if the behavior of &(p 2 ) is given by 
1 — cp 2K + 0(p T ) with t > 4k. Moreover, this quantity goes to 1 at large momenta and its absolute value is clearly 
bounded from above by some positive constant M if a{p 2 ) G [0, 1]. Hence, we have 

?S - C^ 1 -^ 1 + M72 ^°°) = ^M + ^,oc) + M/ 2 ( P 2 ,^) . (148) 

For 1/2 > k we can also use the upper bound in Eq. (B33) and write 

aip2) <(^ + m)i 2 (p 2 ,oo), (149) 



9 2 N C ~ V c 

where M" is a positive constant. Thus, we have again an IR singularity proportional to — ln(p 2 ) if 2?(0) is not zero. 
We conclude that a{p 2 ) can be finite solely if 2?(0) = 0. 

Let us remark that the only hypothesis considered in this case is the IR expansion a(p 2 ) = 1 — cp 2K + 0(p T ) with 
1 > 2k and r > An. Also note that the 2d lattice data [56] show for the ghost propagator an IR behavior in good 
agreement with the so-called scaling solution [14-16] that predicts n — 0.2. Thus, the condition 1 > 2k is verified in 
both cases. One can also note that, by considering in Eq. (139) the full ghost-gluon vertex T\ l/ (p, q), instead of the 
tree- level one 5\ v , the above results still applies for an IR-finite vertex. This hypothesis is usually adopted in DSE 
studies of gluon and ghost propagators [12, 20, 25, 31] and it is confirmed by lattice data [98-101]. 



2. Analysis of the Gribov Form-Factor Using a Spectral Representation 

In this section we analyze the DSE (141) in an alternative way, also avoiding the y-max approximation. To this 
end, let us first consider the ^-integral using contour integration. After setting z — e %e we find 

Jo d6 p 2 [1 - a(p 2 )} =lf dZ z 2 (-q + kz) (k - qz) 1 - a [ (-q + kz) (k - qz)z^ ] ' (150) 

where the integral § dz is again taken on the unit circle \z\ = 1. Clearly, besides the poles at q — k/z and at q = kz 
in the first denominator on the r.h.s. of the above equation, one has to consider possible divergences in the function 

/(Z) = 1 - a[(-q + kz){k-qz)z-^] ' (151) 

In particular, if we assume ghost enhancement, i.e. cr(0) = 1, then f(z) is divergent at z = q/k and at z = k/q. Note 
that these divergences are not necessarily poles of the function f(z). Indeed, f(z) could display a branch cut in the 
unit disc or one passing through it. For example, the usual d — 2 DSE scaling solution has G{k 2 ) ~ \/(k 2 y in the 
limit k 2 — > 0, where v is a fractional number. This behavior signals a non-analyticity for Q{k 2 ) at the origin and 
implies a non-analyticity for the function f(z) at z = k/q or at z = q/k. Also, since the ghost is "massless" we should 
expect that the ghost propagator develops a branch cut along the real axis for k 2 < 0. Then, z = q/k or z = k/q 
would correspond to branch points of the function f(z), making quite difficult the evaluation of the contour integral 
in the above expression. 
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In order to overcome this problem, we make the hypothesis that a spectral representation for the ghost propagator 
can be introduced, i.e. we write 27 



G(P > = -2i Tin = / d< 7T^ ' ( 152 ) 

p 2 \-o(p 2 ) J Q t+p 2 

which reproduces the branch cut in Q{k 2 ) for k 2 < (see for example [102]). If we assume ct(oo) = and write 

i r pit) 



«M- ? y dt TTw (153) 

it is clear that the spectral density p(t) must satisfy the normalization condition 

r oo 

1 = / dt p{t) . (154) 
Jo 

Also note that the tree-level ghost propagator Q(p 2 ) = 1/p 2 corresponds to the spectral density p(t) = 2S(t), where 

S(t) is the Dirac delta function. This case will be used below to recover results obtained in the one-loop analysis 
carried on in Sections II A and II B. In the general case, the spectral density p(t) is proportional to the discontinuity 
of the ghost propagator along the cut. 28 

Considering Eqs. (141) and (152) we can write 

f 27T ,„ l-cos 2 (#) i f°° , „ / , z 2 + z 2 -2 /irr , 

Jo [1 -*(<')] = 4 j, Atp{t) f AZ -p qZ 2 + (p 2 + q 2 + t) Z -p q (155) 

i r~,«i** (* a -D a (156) 



/*°° r 

Jo dtp{t) J ' 



4pq Jo ' X ' J z 2 z 2 - (p 2 + q 2 + t) zj (pq) + 1 ' 

where we indicated with z the complex-conjugate of z = e l6 ' . Thus, using the representation (156) we can avoid dealing 
directly with the integral of an unknown function along the branch cut. In exchange, we have in our formulae an 
extra integration of the (also unknown) spectral density p(t). Nevertheless, as we will see below, the above equation 
will allow us to control the p 2 — > limit [at least in the case p{t) > 0]. To this end, let us first note that in the contour 
integral (156) there is a double pole at z — and there are single poles at 



{f + q 2 +t)± J(p 2 + q 2 + tf - 4p 2 q 2 

z ± = ^— VI . (157) 

2pq 

Since p, q, t > we have that p 2 + q 2 +t > 2pq > and one can check that the pole z_ lies within the unit disc while 
z + lies outside of it. Moreover, for p 2 + q 2 + 1 = 2pq (which implies t = and p = q) the two poles coincide and we 
have z± = l. It is also easy to check that the residues, inside the unit circle, for the z-integrand are 

p 2 + q 2 +t 

TZes z=0 = 2~~2 ' ( 158 ) 

pZgZ 

pZqZ 

Then, using the residue theorem, we find 

]„ s 2 [1 - »(s 2 )] 2 1„ " p 2 q 2 



Since we are working in the d = 2 case, the theory should be UV finite and we do not need to consider rcnormalization factors here. 
Note that, if Q(k 2 ) has a branch cut along a curve C in the complex plane and if it goes to zero sufficiently fast at infinity, using 
Cauchy's theorem we could write down an integral relation similar to Eq. (152) with the variable t running over the curve — C, with 
z 6 — C <t=> — z S C. Also, possible poles can be included by adding <5-functions to the spectral density p(t) or, equivalently, by pulling 
the pole terms out of the spectral integral. 
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and we can write the ghost DSE (141) as 

5 2 7V C 7 8tt vy ; 7 FK) p 2 q 2 

+ x + t - ^Jt 2 + 2t (p 2 + x) + (p 2 - xf 



Note that, for p(t) = 25(t) and using 



VKp + v) 2 } [(p-q) 2 } = 



p 2 x 



p 2 - q 2 if p 2 > q 2 

q 2 - p 2 if q 2 > p 2 



(161) 
(162) 



(163) 



we find from Eq. (160) the one-loop result (11). Also note that, by Taylor expanding the integrand at p 2 = 0, one 
finds 



*(0) 
g 2 N c 



1 f°° j nf . [°° j pit) 1 °° , Z> x 
= — / dxX>(x) / dt-^- = — / dx— — V — L r —. r , 
8tt i w Jo t + x 8tt 7o x 1-ct x 



(164) 



where we used the definition (152). As shown above [see Eq. (144)], this result can also be obtained immediately 
by setting p 2 = in Eq. (141). However, as already pointed out below Eq. (144) and in Ref. [89], one should not 
exchange the q integration and the p 2 — > limit. Therefore, in order to properly evaluate a(p 2 ) for small momenta 
p 2 , we write Eq. (162) as 



a(0) 



= lim 



g 2 N c P *->o J 16tt 



dx f°° p 2 +x + t - Jt 2 + 2t(p 2 + x) + (p 2 - xf 
V{x) / dtp(t) i 



p 2 x 



r°° dx r°° 

lim / — —V(x) / At pit) 

p*->0 J„2 167T V ' J ^ ' 



p 2 +x + t - \jt 2 + 2t(p 2 +x) + [p 2 -xf 



p 2 x 



(165) 



The first integral can be estimated using the the trapezoidal rule. We then obtain 

r°° p 2 + x + t - \jt 2 + 2t (p 2 + x) + (p 2 - xf 



= lim 

p 2 ^0 



f 


dx 


Jo 


167T 


P 2 


f- QO 


32tt 


Jo 


1 
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Jo 



p 2 x 



dt p(t) 



, A 2p 2 +t - yjt 2 + Up 2 2X>(0) 

p* + ~ P T Tt 



v( P 2 ) 



2p 2 + t - yjt 2 + Up 2 



+ lim 



2>(0) 



p^o 16tt[1 - crip 2 )} 



(166) 
(167) 

(168) 



where we used again Eq. (152). For the second integral we define 

j 2 + x + t - \jt 2 + 2t ip 2 + x) + ip 2 - xf 



G(x,p ) - 



16?r J 



dt pit) 



and find 



j dx Q ( X,P ' = ln(x) Gix,p 2 ) \°° 2 - [°° dx ln(x) G'ix,p 2 ) 

2p 2 + t - ^Jt 2 + Up 2 



ln(p : 



M V(P 2 ) ' 
16tt 



dt p[t) 



dx ln(x) g'ix,p 2 ^ 



(169) 

(170) 
(171) 
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where ' refers to the derivative w.r.t. the x variable and we used the fact that T>(x) goes to zero at large x. Note that, 
in the one-loop case p(t) — 2S(t), we have Q(x,p 2 ) = V(x)/(8ir) and Eq. (171) becomes 



dx V{x) 
,2 8ir x 



= - Hp 2 ) 



v(p 2 ) 

8tt 



Jp 2 



dx 
8^ 



ln(a:)I>'(a:) , 



in agreement with Eqs. (28) and (30). 

By collecting the above results we can therefore write 



^ = J_ limj P W 
g 2 N c 16tt p 2 ^o | 1 _ a ( p 2) 



V(P 2 ) 



Hp 2 ) 



dtp(t) 



2p 2 + t - ^Jt 2 + Up 2 



f 

Jv 2 



dx \n(x)G'(x,p 2 ) } . 



We can now verify that the last integral in the above expression is finite. Indeed, we have 



G'(x, P 2 ) 



+ 



16tt 
16tt 



dt p{t) 

Jo 



p 2 +x + t - 



2t (p 2 + x) + (p 2 - xf 



x + t — p 2 



Then, for large x we find 
G'{x lP 2 ) 



8tt 



f 

Jo 



dt p{t) 



p2 p 2 ^t 2 + 2t{p 2 + x) + (p 2 -x) 
V{x) f 



8irx 2 



dt p(t) t 



8tt 



dt p(t) t 



where we used the normalization condition (154), while for small x we have 



G'(x,P 2 ) 



xV'jx) 
8tt 



dt 



pit) , V(x) 



t+p 2 



8tt 



dtM; = ^[xV>(x)+Vix)} 



t + p A 



8tt 



(172) 



(173) 



(174) 



(175) 



(176) 



where we used the definition (153). Thus, the integral J °° dx ln(x) G'ix,p 2 ) has no IR and UV singularities if T>(x) 
and T>'{x) goes to zero sufficiently fast when x goes to zero and to infinity. At the same time we need the integral 



Jo 



dtpit)t 



to be finite. We can also check that the integral 



f 

Jo 



dtp(t) 



2p 2 + 1 - yjt 2 + Up 2 



(177) 



(178) 



is finite and nonzero if p(t) is non-negative. 29 To this end let us first note that the numerator 2p 2 + t— \/t 2 + Atp 2 

(182) 



is non- negative since 2p 2 + t > \Jt 2 + Atp 2 when t,p 2 > 0. Moreover, if we define 

$(t,p 2 ) 



2 P 2 + t - yjt 2 + Up 2 



29 Our results cannot be easily extended to the general case of a spectral density p(t) that is negative for some values of t. However, they 
apply if one can explicitly verify that the integral in Eq. (178) is indeed finite and nonzero. Also note that we cannot simply consider 
the limit p — > of the integral (178), since the factor multiplying p(t) vanishes in this limit and we might erroneously conclude that the 
above expression is equal to zero. Indeed, already in the tree-level case, i.e. for p(t) = 28(t), one finds that the integral (178) is non-zero 
and equal to 2. Finally, one could also formally expand the integrand in powers of p 2 leading to 

2p 2 +t - ^t 2 + Up 2 _ ™ { 1/2 \ f4p 



where ( a ) = — ^ Q 7- ^ r is the usual binomial coefficient. However, this series is not converging for all values of t 6 [0, +00). 

\b J r(6 + l)r(o-6 + l) 

Moreover, for n = 2 we have a term proportional to 

r dt (ISO 

Jo t 

and this integral is divergent. Indeed, by (formally) setting p 2 = in Eq. (152) we find 

'dt^=oo. (181) 



■/() 
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it is clear that 2 = $(0,p 2 ) > &(t,p 2 ) > 0, since the quantity t — \/t 2 + Atp 2 is negative for t > 0. This implies 

f-OO />oc 

/ dtp(t)$(t,p 2 ) < 2 / cltp(t) = 2, (183) 
Jo Jo 

where we used again the normalization condition (154). At the same time we can write 



Jo 7o * + P 2 P 2 2 7 



d^, (184) 

where we use the fact that the function (t + p 2 )Q(t,p 2 )/p 2 is positive and gets its minimum value, equal to 3/2, for 
t/p 2 = 1/2. Then, using the definition (153) and the condition cr(p 2 ) > 0, we can write 

r°° 33 

I dt ^*^>w^m\-i- (185) 

From the above results we conclude that in Eq. (173) we have two possible IR singularities, i.e. the term 2?(0)/[l — 
<j{p 2 )] 1 if er(0) = 1, and the term proportional to —V(p 2 ) ln(p 2 ). In both cases we need to impose the condition 
T>(0) = in order to avoid the singularity. Thus, we find again that a massive gluon propagator in the d = 2 case is 
not compatible with the restriction of the functional integration to the first Gribov region. 



B. The 3d Case 

In the 3d case we expect no UV divergences when using dimensional regularization 30 and the DSE for the Gribov 
form-factor is simply 



a{p 2 ) 1 f d 3 q „, 2 , p n l 1 

= p^ J j2^P^M avT> (9 )P^(q)^Y^{ 



(186) 



g 2 N c p 2 J (2tt)3^ " v "v* > ^ v ^> s>\-<j{s 2 ) 

{ £r V{q2) I dn3 ^[l-^s 2 )] 
where we used the tree-level ghost-gluon vertex T\ v (p, q) = 5\ v and s 2 = p 2 + q 2 — 2pq cos(<^i). We can now work as 
in Section IV A 1 and use the results of Appendix B. In this way we obtain the upper bounds 



o 



g (P 2 ) < h{p 2 ,oo) nRR . 
-JN C ^ 1-7(0) ' (188) 



^ <(^+M)W,oo) t (189) 



if a(p 2 ) < cr(0) < 1, and 

9 2 N C V c 

if a(p 2 ) < ct(0) = 1 with <r(p 2 ) w 1 — cp 2K + 0(p T ) for small p 2 . In the latter case we also need the conditions 1 > k 
and r > 4k. As we saw in Eq. (B35), under simple assumptions for the gluon propagator T>(q 2 ), the integral ^{p 2 , oo) 
is finite in the IR limit p — > 0. Thus, in both cases the upper bound of a(p 2 ) is also finite 31 and, in order to have 
a finite value for a(0) in the 3d case we do not need to set T>(0) = 0. This result also applies when an IR-finitc 
ghost-gluon vertex is included in the ghost DSE (186). Let us also note that the scaling solution predicts in the 3d 
case [14-16] a value k w 0.4 for which the condition 1 > n > is satisfied. 



As shown in Section IIIB above at one-loop level, the evaluation of the ghost propagator in 3d usually involves Gamma functions with 
half-integer arguments, which do not generate infinities. Indeed, for nonnegative values of n with n integer, one has [97] T (n + 1/2) = 
y^2- n (2ra - 1)!! and T (-n + 1/2) = (-2)™ ^/ (2n - 1)!! , where n\\ denotes the double factorial. 
Using the fact that <r(p 2 ) is nonnegative and Eqs. (B31) and (B36), the lower bound 

^/ 3 (p 2 ,oo) < 7(p 2 ,l,3,oo) < (190) 
4 g 2 N c 

clearly applies. However, since c(p 2 ) is finite, this bound does not add any relevant information to our analysis. Note that this 

observation can be made also for the Ad case described in Section IV C. 
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C. The 4d Case 



In Ad the DSE for a{p 2 ) is given by (see for example [94]) 

«f) - i - a + 2. jT d, ^ vch I «. p'f^ff)] • (191) 

where Z 3 and -Zi are the renormalization constants for the ghost propagator and the ghost-gluon vertex respectively 32 
and s 2 = p 2 + q 2 — 2pq cos(^i). In order to eliminate these constants from the expression for <r(p 2 ) we can subtract 33 
the same equation for some fixed value p 2 = /i 2 and set Z\ = 1, using the non-renormalization of the ghost-gluon 
vertex in Landau gauge [103]. This gives 



a(p 2 ) a{ P 2 ) [°° g 3 y l-cos 2 ^) 

?^-^: + io dq wr v{q) J dni s 2 a - a(s 2 )] 



-r«j&™j^*r£%> (i92) 



with t 2 = fi 2 + q 2 — 2fiq cos(<^i). Clearly, considering T>(q 2 ) <~ 1/g 2 at large momenta, the UV (logarithmic) 
divergence of the two integrals cancels out. In order to derive upper bounds for the above expression without spoiling 
the cancellation of UV divergences and since we are interested in the IR limit p — ► 0, we can consider a momentum 
scale t sufficiently large and write 

a(p 2 ) a(^) f* g 3 f l-cos 2 (<M 

+ f I ^^tk - f d « (4^ / ^^0M-y <™ 

For small momenta p only the first integral on the r.h.s. of the above equation can produce an IR singularity. Following 
the analysis presented in the 3d case above we can then write 

, <? 3 ^/ 2^ /* ,^ l-COS 2 (0i) / 4 (P 2 ,^) 



if <y(p 2 ) < cr(0) < 1, and 



/ ^ ^^^l < f — + M ) 7 4 ( P V) , (195) 



(2tt) 4 ^ ' J s 2 [1 - ^(s 2 ); 

if u(p 2 ) = 1 — q? 2K + 0(p T ) with t > 4k and 3/2 > n. Again, thanks to the result (B35), both upper bounds are finite 
in the IR limit p -> also for X>(0) > 0. 

An alternative proof can be given by working directly with Eq. (191) and using dimensional regularization, i.e. 
considering a dimension d = 4 — e. In this case we can write 

a(p 2 ) = 1 - Z 3 + Z x <r d (p 2 ) (196) 

with 34 

r iq ^ {f) i iai (197) 



J„ *(2ir)< J j2[l-^(» 2 )] 



Note that we are again considering a tree-level ghost-gluon vertex r A „(p, g) = 8\„. 
Again, this corresponds to a MOM scheme with the condition G(p, 2 ) = 1//J. 2 - 

Of course, with d = 4 — e and e > 0, the integral in Eq. (197) is no longer dimcnsionless. To keep the dimensionality correct we should, 
as always, scale out a dimensional factor m 4 ~ d where m is a mass scale, which could then be combined with the coupling constant g 2 , 
making cr(p 2 ) dimcnsionless also for e > 0. This is important when evaluating the e expansion in order to single out UV divergencies. 
Since here we are mainly interested in the IR behavior of <r(p 2 ), we do not keep track explicitly of all the terms depending on e and we 
simply consider the coupling g 2 dimensionful. 
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Then, if we can show that no IR singularities occur for 2 < d < 4, the UV infinity that appears for d — > 4 is taken 
care of by the renormalization factors. In order to show that <Jd{p 2 ) is IR finite we can work as done above and write 



adip2) < ( ¥L + M ) , (199) 



or 

9 2 K 

depending on the value of Cd(0). In the latter case we considered again the IR expansion (Jd{p 2 ) ~ 1 — cp 2K + 0(p T ) 
and the conditions r > An and 3/2 > k . We conclude that also in the Ad case, tr(0) is finite if 2?(0) is also finite (but 
not necessarily null). 

Let us remark that the IR exponent usually obtained in the scaling solution [14-16] is n « 0.6 in the Ad 
case, i.e. the condition k < 3/2 is satisfied. Also note that when (d — l)/2 < k the hypergeometric function 
2F1 (1 + k, 1 + k — d/2; 1 + d/2; z) is not convergent at z = 1 and we cannot use the above proofs in order to derive 
properties of the Gribov form-factor. However, these large values of k imply for the ghost propagator G(p 2 ) a very 
strong IR enhancement with a behavior at least as singular as 1/fc 5 in Ad and at least as singular as 1/fc 4 for d = 3. 



V. CONCLUSION 

Summarizing, in this manuscript we have considered general properties of the Landau-gauge Gribov ghost form- 
factor a(p 2 ) for SU(_/V C ) Euclidean Yang-Mills theories in d > 2 space-time dimensions. This form-factor is in a 
one-to-one correspondence with the ghost propagator G{p 2 ) via Eq. (2). Also, a(p 2 ) is bounded by 1 if the no-pole 
condition (3) is imposed, i.e. if one restricts the functional integration to the first Gribov region ft. The main result of 
this work is an exact proof of qualitatively different behavior of a(p 2 ) for d — 3, 4 with respect to d = 2. In particular, 
for d = 2, the gluon propagator T>{p 2 ) needs to vanish at zero momentum in order to avoid in u{p 2 ) an IR singularity 
proportional to — 2?(0) ln(p 2 ). On the contrary, for d = 3 and 4, an IR-finite ghost form-factor a(p 2 ) is obtained also 
when T>(0) > 0. These results were proven, in Section II, using perturbation theory at one loop and, in Section IV, 
by considering the DSE for the ghost propagator. Let us stress again that in DSE studies of correlation functions in 
minimal Landau gauge, besides using the no-pole condition, a specific boundary condition is usually imposed on the 
Gribov ghost form-factor at zero momentum. Here, instead, we have tried to prove general properties of the Gribov 
ghost form-factor a(p 2 ) when the restriction to the first Gribov horizon is considered. 

At the same time, in Section III, we have presented closed analytic expressions for the Gribov form-factor a(p 2 ) 
at one loop, considering for the gluon propagator linear combinations of Yukawa-like propagators (with real and/or 
complex-conjugate poles). These functional forms, briefly described in Eqs. (53)-(55), were recently used to fit lattice 
data of the gluon propagator in the SU(2) case [59, 60]. The expressions obtained for <r(p 2 ) confirm the results 
presented in Section II. These expressions also show that, for the ghost propagator G(p 2 ), there is a one-parameter 
family of behaviors [21, 34] labelled by the coupling constant g 2 , when it is considered as a free parameter. The 
no-pole condition <r(0) < 1 then implies g 2 < g 2 , where g 2 is a "critical" value. For g 2 smaller than g 2 one has 
<r(0) < 1 and the ghost propagator is a massive one. On the contrary, at the "critical" value g 2 , i.e. for a(0) = 1, one 
finds an IR-enhanced ghost propagator. As stressed in the Introduction, the physical value of the coupling is expected 
to select the actual value of c(0). Present results [21, 34] give ct(0) < 1 in the four-dimensional SU(3) case. 

Our findings imply that a massive gluon propagator cannot be obtained in the two-dimensional case, in disagree- 
ment with some of the results presented in Ref. [16] (see their Table 2). 35 A possible massive behavior for the gluon 
propagator in the 2d case was also explicitly conjectured in Ref. [65] as a Gribov-copy effect. However, since our 2d 



A massive solution in the 2d case was also obtained in Ref. [14]. On the other hand, the author of [14] stressed that a full understanding 
of the 2d case would require a more detailed investigation. 
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result is valid for any Gribov copy inside the first Gribov region, we have shown that, at least for the 2d gluon propa- 
gator V{p 2 ) in the minimal Landau gauge, Gribov-copy effects do not alter our conclusion for the value of the gluon 
propagator at zero momentum, i.e. 2?(0) must vanish. This observation also represents an explicit counterexample to 
the identification of the one-parameter family of solutions for the gluon and ghost DSEs [22, 31] with different Gribov 
copies, as suggested in [65, 67, 68]. 
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Appendix A: Angular Integration in the d-Dimensional Case 

Following Appendix B in Ref. [15] one can easily perform, using hypcrsphcrical coordinates, the angular integrations 
necessary for our calculations. To this end let us recall that, in d dimensions, one has the following relations between 
Cartesian coordinates Xi (with i = 1, 2, . . . , d) and hyperspherical coordinates r, cfij (with j — 1,2, ... ,d— 1): 

xi — r cos(</>i) , 

x 2 = r sin(0i) cos(</> 2 ) , 

x 3 = r sin(0i) sin(<p 2 ) cos(^ 3 ) , 

x d -i = r sin(0i) sin(0 2 ) • • ■ sin(^ d _ 2 ) cos(0 d _i) , 

x d = r sin(0i) sin(</> 2 ) . . . sin(</> d _ 2 ) sin(<^_i) . (Al) 

The hyperspherical coordinates take values, respectively, r <G [0, oo), <j>i <E [0, w] for i = 1, 2, . . . , d—2 and <fid-i € [0, 2ir). 
At the same time, the volume measure is given by 

dy = r d " 1 drdO d (A2) 

with 

dfl d = sin d - 2 (<Ai) sm d - 3 (<j> 2 ) . . . sin(0 d _ 2 ) d& d</> 2 . . . d0 d _ 2 d<^_i . (A3) 

Note that, in the usual three-dimensional case, this notation correspond to x\ = z, x 2 — x and x 3 = y. 
Here we want to evaluate the integral 



/ 



dfi d /(0i) • (A4) 

If wc indicate with Q d the well-known result 



2w d / 2 

r(|) 

where T (x) is the Gamma function, then we have 



= J dn d = —^, (A5) 



/ /(&) dn d = , d n 2 d f S in d - 2 (^)/(^i) dfa 

J J sin (0i)d0i Jo 



(A6) 



(•7T . 

Jo sir 

where all other angular integrations have already been evaluated. The integral in the denominator can be written as 

'd-1 V 



P sin^OOd^ = f (l-z 2 f- 3)/2 dz = f 1 t-W (l-t) (d - 3)/2 dt = B 
Jo J -1 Jo 



2 ' 2 



(A7) 



31 



where B (a, b) = T(a) T(b)/T(a + b) is the Beta function. 

In our calculations we consider two different functions / (0i), i.e. 

f (<l>i) = 1- cos 2 (00 (A8) 

and 

f(A\ 1 - COS 2 (01 ) 

[p 2 + q 2 - 2pq cos(0i)] 

In the first case the integration gives 

f [l-cos 2 ^)] dn d = n sm d - 2 (^) [l-cos 2 ^)] d0i 

J « ( — , 2J Jo 

and the integral in the numerator is 

f (1 - z 2 ) (d - 1)/2 dz = f t-V" (1 - i)^" 1 )/ 2 dt = B (^±1 1 
Collecting these results we find 



we have 



(A10) 



(All) 



/ [1 - cos 2 (^)] dO d = O d Z/ 2 = fi d — — , (A12) 

J B \ — >2) a 

where we used xT(x) = T(x + 1). In the second case, i.e. when considering the integral 

/ -7^1(00 r dQd ' (A13) 



v #1 • (A14) 



Jo [p 2 + <? 2 - 2 Pg cOs(0l)] 
We can now use the result (see for example formula 3.665.2 in [97]) 

I [i + fl2±2fl cos W r dd = 1/2) 2Fl ^ u ~" + 1/2; " + 1/2; fl2) ' (A15) 

which is valid for \a\ < 1 and Re(/i) > 0. Here 2-F1 (a, c; 2) is the Gauss hypergeometric function (see Appendix B). 
Then, we find 

/ r 2^ V'f^ ^^ dnd = \ n \ B (^>k) 2F 1 (,,,-d/2;l + d/2-p 2 /q 2 ) (A16) 

J [p 2 + q 2 - 2pq cos(0i)] B , 3) 9 V 2 2 / 



^ ^ 2*1 (f, * - d/2; 1 + d/2-p 2 /q 2 ) (A17) 



if p 2 < q 2 and 

/ f 2^ ^""f^ = ~j ~r~ - - d / 2 ; 1 + d / 2 ; ? 2 /p 2 ) (Ms) 

7 [P + 1 — 2pqcos{4>\)\ p 2 a v y 

if g 2 < p 2 . 
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Appendix B: Properties of the Gauss Hypergeometric Function 



Let us recall that the Gauss hypergeometric function 2 .Fi(a, b; c; z) is defined [97] for |z| < 1 by the series 

j?( h ^ («)» (b)n z n . , ab a(a + 1)6(6 + 1) z 2 mn 

where 

_ r(a + n) 

(a) - - ( } 

is a so-called Pochhammcr symbol. This series is converging for c ^ 0, — 1, —2, .... It is also converging for \z\ = 1 if 
5i(c — a — b) > 0, where 3? indicates the real part. When this condition is satisfied one can use, for z — 1, the result 
(see formula 9.122.1 in Rcf. [97]) 

■^^■'- ^-l (B3) 

In Eqs. (A17) and (A18) the hypergeometric function appears with c = l + d/2. Therefore, the corresponding series 
is converging inside the unit circle for any dimension d > 0. At the same time we have c — a — b = d + 1 — 2v and 
2 F\(v, v — d/2; 1 + d/2; z) is finite on the unit circle for (d + l)/2 > v. Also, using the relation (see equation 9.131.1 
in Rcf. [97]) 

2-Fi (a, b-c-z) = (1 - z) c - a - b 2 iq (c - a, c - 6; c; z) (B4) 

we can write 

2 F X {v, v - d/2; 1 + d/2; z) = (1 - z ) d+1 - 2u 2 *i (1 + d/2 - v, 1 + d - v; 1 + d/2; z) (B5) 

and for (d+l)/2 — f > we have that l+d/2 — v, 1+d—v > 0, i.e. the hypergeometric function 2 *i (^, v—d/2; l+d/2; z) 
is clearly positive for z € [0,1]. 

In the case v = 1 the above results simplify. In particular, we have convergence of the series on the unit circle for 
any dimension d > 1. Then Eq. (B3) yields 

Fni rf/9-i i rf/o-n - r(i + rf/2)r(d-i) d/2 

2Fl (1 ' 1 ~ rf/2 ' 1 + rf/2 ' 1} - r(d)r(d/2) - J— i ' (B6) 

which is positive. Note that for d = 2 the above result gives 2 F\ (1, 0; 2; 1) = 1. One can actually check that, when 
d = 2, the function 2 *i (1, 1 — d/2; 1 + d/2; z) = 2 F\ (1, 0; 2; z) is simply equal to 1 for any value \z\ < 1. Indeed, 
from the series representation (Bl) it is obvious that 2 Fi(a, b; c; z) — 1 for b — (and/or for a — 0). The same result 
can be obtained by considering Eq. (B4) and the relation (see equation 9.121.1 in Rcf. [97]) 

2 Fi (n,c;c; s) = (1 -*)-", (B7) 

yielding 

2 F 1 (a, 0; c; z) = (1 - z) c ~ a 2 iq(c - a, c; c; z) = (1 - zf- a (1 - z) a ~ c = 1 . (B8) 

The hypergeometric function 2 F\ (1, 1 — d/2; 1 + d/2; z) is actually very simple for any even dimension d = 2n, with 
n > 1. Indeed, in this case we are considering the function 2 Fi (1, 1 — n; 1 + n; z) and the coefficient 6 is either zero 
or negative. This implies that the series (Bl) is actually a polynomial in z. For example, for d = 4 we have the very 
simple expression 

2^(1,-1; 3; z) = 1 - | . (B9) 
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Note that it is also possible to find a closed form for the hypergeometric function 2-F1 (1, 1 — d/2; 1 + d/2; z) in the 3d 
case. Indeed, by considering the relations (see equations 9.137.8 and 9.137.1 in Rcf. [97]) 

= 2 F 1 (a,b;c;z)c + 2 F 1 {a + 1, b; c + 1; z) (6 - c) - 2*1(01 + 1, b + 1; c + 1; z) b (1 - z) (BIO) 

= 2-F1 (a, b; c; z) c [c — 1 — (2c — a — b — 1) z] + 2*1 (a, 6;c+ 1; z) (c - a) (c - b) z 

+ 2 F 1 (a,b;c-l;z)c{c-l) (z-1) (Bll) 

and (see equation 9.131.1 in [97]) 



we can write 



2 Fi(a,&;c;z) = (1-z) 2*1 ( 6, c - a; c; — - ) (B12) 



2F1 (1, -1/2; 5/2; z)=l 2 F 1 (0,-1/2; 3/2; z) + 2 F X (1, 1/2; 5/2; z) (B13) 

= i + ^^7^ [2^1 (1,1/2; 1/2; z) - 2 Fi (1,1/2; 3/2; z)] (B14) 
= I + [1 ~ (!-^) 2^1 (1,V2; 3/2; z)] (B15) 



= 3 3(l-z) 
4 8z 

= 3 3(l-z) 
4 8z 



1 - JY~z 2 F x (\/2, 1/2; 3/2; 



z- 1 



1-i^arcsinh 



\/z VV 1 - z 



(B16) 
(B17) 



where we have also used Eq. (B7), the relations 2*1(0, b; c; z) = and (see 9.121.27 in [97]) 



2*1(1/2, 1/2; 3/2; -z 2 ) = ™5*£ . (B18) 

z 

From the expression (B17) it is easy to check that 2 *i(l, -1/2; 5/2; 1) = 3/4 and that 2 *i(l, -1/2; 5/2; 0) = 1, as 
expected. 

Using the above series (Bl) one can verify that the derivative of 2*1 (a, b; c; z), with respect to the variable z, is 
given by 

— 2 * 1 (a,&;c;z) = — 2 *i(a + 1, b + 1; c + 1; z) . (B19) 
oz c 

Thus, in the case of interest for us, we have 

^ 2 F x {y,v - d/2;\ + d/2;z) = " 2*i(^ + 1, ^ + 1 - d/2; 2 + d/2; z) . (B20) 

These results can be written as 

^ 2*1 (1/, 1/ - d/2; 1 + d/2; z) = - ^ ~ ^ (1 - zf~ 2 " 2*1 (1 + d/2 - v, 1 + d - v; 2 + d/2; z) , (B21) 

where we made use of Eq. (B4). When the hypergeometric function 2 Fi(i/,v — d/2; 1 + d/2;z) is finite in the unit 
circle, i.e. for (d+ l)/2 > v, it is clear that this derivative is positive for v > d/2 and 2 Fi{v, v— d/2; 1 + d/2; z) attains 
its maximum value at z = 1, equal to 
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and its minimum value, equal to 1, at z = 0. On the contrary, the same derivative is negative when v < d/2. 
In this case 2 Fi(f,v — d/2; 1 + d/2;z) is largest at z = 0, with 2 F\(y,v — d/2;l + d/2;0) = 1, and smallest at 
z = 1 with a value given by Eq. (B22). Finally, for v = d/2 the derivative is null and the hypergeometric function 
2 F\(v, v — d/2; 1 + d/2; z) is equal to 1 for all values of z e [0, 1]. Thus, for (d + l)/2 > v we can always write 
M' < 2 F\(y,v — d/2; 1 + d/2;z) < M" , for some positive constants M' and M" and with z taking values in the 
interval [0, 1]. 

For v = 1 these results again simplify, yielding 

^- 2 F 1 (l,l-d/2;l + d/2;z) = 2 iq(2, 2 - d/2; 2 + d/2; z) . (B23) 

As expected, for d = 2 this derivative is zero since 2*1 (1,0; 2; z) = 1. The result above also simplifies for d = 4, for 
which we find on the right-hand side of Eq. (B23) the value — 2*1 (2, 0; 4; z)/3 = —1/3, as already known from Eq. 
(B9). At the same time, for d = 3 we can write 

^ 2*1(1,-1/2; 5/2; z) = -\ 2 *i(2, 1/2; 7/2; z) = -\ (l + ^ + ^ (B24) 

and the derivative is clearly negative for any value z > 0. The same result can actually be proven for any dimension d 
larger than 2. Indeed, the hypergeometric function 2-F\(2, 2 — d/2; 2 + d/2; z) is finite in the unit circle for d > 2 and 
using Eq. (B4) we can easily verify that it is also positive for z € [0, 1]. Thus, the derivative in Eq. (B23) is negative 
for d > 2 (and z e [0, 1]). As a consequence, under the same hypotheses, we have that the hypergeometric function 
2*1 (1, 1 — d/2; 1 + d/2; z) has its maximum value, equal to 1, at z = 0, and its minimum value, equal to d/ (2(d — 1)), 
at z = 1 [see Eq. (B6)]. 

Using the definition (B2) of the Pochhammer symbol (a) n it is also easy to verify that 

(a) n = a (a + 1) (a + 2) . . . {a + n - 1) (B25) 

which implies 

n{a) n = a [(a + l)„ - (o)„] . (B26) 
Then, using Eq. (Bl), one can prove the relation [104] 

z-2* 1 i( a ,6;c;z)^^^— n- (B27) 



n=0 



(q + 1)„ (&)„ Z™ _ y, (q)„ (&)„ Z^_ 
_ n— n— 



(B28) 



= a [ 2 Fi(a+l,6;c;z) - 2^1(0, b; c; z) ] . (B29) 

In particular, we can write 

z-^- 2 Fi(l,l-d/2;l + d/2;z) + 2 *i(l, 1 - d/2; 1 + d/2; z) = 2 *i(2, 1 - d/2; 1 + d/2; z) . (B30) 
oz 

Note that the r.h.s. in the above relation is finite in the unit circle for d > 2 and, using again Eq. (B4), we can verify 
that it is also positive for z e [0, 1]. 

The above results, together with the Eqs. (A17) and (A18), allow to write a lower and an upper bound for the 
integral 

I(p 2 ,v,d,l) = f dq^Vtf ) ( An d \- C f {<t,l) (B31) 

Jo ( 27 J [p 2 + q 2 -2pqcos((j> 1 )} 
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Indeed, after considering the angular integration, we have (for t > p) 



i{ P 2 ,v,d,t) 



n rl d-i 



(2TT) d d 







P d l l^ 1 1 ^r L ^ - d/2; 1 + d/2; q 2 / P 2 ) 



+ / dqq d - 3 V(q 2 ) 2 F 1 (v,v-d/2;l + d/2;p 2 /q 2 ) 



Then, for 1 + d — 2v > we obtain 



M'I d (p 2 ,£) < I( P 2 ,v,d,£) < M"I d (p 2 ,£) 



with 36 



d- 1 



(27r) d d 

Also, in the limit p — > 0, we find (for d > 1) 
lim/ d (pV) = 



lim £!« 

p^o 2 



dg g d - 3 2?(g 2 ) 



/' 

Jo 



dqq d ~ 3 T>(q 2 ) 



(B32) 



(B33) 



(B34) 



(B35) 



where we used the trapezoidal rule. Clearly, for T>(0) > 0, the first term is IR-finite if d > 2 while the second term is 
finite for d > 2. Finally, note that for v = 1 the condition 1 + d — 2v > simplifies to d > 1. At the same time, for 
> 2 the inequalities (B33) become 

d 



2(d-l) 

and for d — 2 we have I(p 2 ,v, 2, £) = ^(p 2 , (■)■ 



I d {p z ,i) < I{p 2 ,v,d,t) < I d {p 2 ,£) 



(B36) 



Appendix C: Hypotheses on the 2d Gluon Propagator 

In Section II B we have proven, using two different approaches, that in the 2d case one needs to set T>(0) = in 
order to have cr(p 2 ) < +oo. The assumptions made for the gluon propagator were rather general. Indeed, for the first 
proof one needs, for small momenta p 2 , an expansion of the gluon propagator of the type V(p 2 ) w 2?(0) + B p 2r > + Cp 2( - , 
with £ > r/ > and T>(0), B and C finite. At the same time, for large momenta p 2 , we required 

lim V(p 2 ) = lim = . (CI) 

p 2 — >oo p 2 — >oo p 

Let us recall that we are indicating with D{p 2 ) a primitive of V{p 2 ) and that V{p 2 ) is the first derivative with respect 
to the variable p 2 . In Section II B above we considered for T>(p 2 ) a large p 2 behavior of the type 1/p 2 - However, it is 
clear that a weaker condition can also be used. Indeed, the behavior T>(p 2 ) <~ 1 /p 2e with 1 > e > also allows to satisfy 
the above conditions. In order to check this, one should recall that T>(p 2 ) <~ l/p 2e implies 37 £>{p 2 ) ~ p 2 ~ 2t + constant 
and these two asymptotic behaviors yield the limits in Eq. (CI). Under the same hypothesis we can also verify that 
the integral [see Eq. (26)] 

ix l^ + M)- M < D ^-^\-^ x " + M ^ (C2) 



Jo 



Note that, for a gluon propagator V(q 2 ) with a behavior 1/q 2 at large momenta, the second integral in Eq. (B34) is UV divergent if 
d > 4 and t = oo. 

Let us recall that, while this (Abelian theorem) is a correct statement, the converse (also called Tauberian theorem), i.e. D(p 2 ) ~ p 2 ~ 2e 
implies T>(p 2 ) ~ 1/p 2e , is not always true (see for example Ref. [105] and Section 7.3 in Ref. [106]). This is why the so-called de 
l'Hopital's rule does not always apply (see also footnote 39). 
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is finite. To this end, let us first consider the integral 

/ 

Jo 



" dx in W + M) vMiV^-VjO)] X ^ + M)V'( X ) 



with < p and p small. In the limit x — > the integrand behaves as 

ln(M) [MC (n-£) x^- 1 - Brtx 71 - 1 } 



(C3) 



(C4) 



and no singularity arises in the integration from x = to x = p 2 (if £ > 77 > 0). At the same time, for large x the 
same integrand behaves as 



ln(x) 



#P(0) 



~ V'{x) 



Thus, for sufficiently large £ 2 , the integral 



dx In (a" + M) 



rjM [D(x) - £>(0)] - x (x^ + M) P'(x) 



can be approximated by 



/oo 
dx ln(x) 
, 2 



r)MV(0) 
x 1+r i 



V'{x) 



After integrating by parts 38 we then find 



I a (£ 2 ) = -ln(x) 



MX»(0) 



xi 



V(x) 





00 

+ 


f-OC 

/ dx 









MV(0) V(x) 



(C5) 



(C6) 



(C7) 



(C8) 



which is clearly finite under the assumptions made for the gluon propagator V(p 2 ). Finally, the remaining term 

Jl dx In W + M) nM[V(x)-V m ^x(x^ + M)V' { x) ^ 



v 1 x l+ n 

with p<Cl, can be easily bounded by making the assumption that neither V{x) nor V'(x) display a singularity for 
x e [p 2 J 2 ]. Thus, we can conclude that the integral (C2) is indeed finite. We further notice that I a (£ 2 ) is null in the 
limit £ 2 — > 00. Going back to the first proof in Section II B, one can verify that the above conditions allow to show 
that the Gribov form-factor a(p 2 ) goes to zero for p 2 going to infinity and that the term —T>(0) ln(p 2 ) is the only 
singularity of a(p 2 ) in the IR limit p 2 — > 0. 

The situation is, of course, very similar in the second proof. In this case we considered a finite value for X>(0) and 
the limit 



lim [p 2 \n(p 2 ) - p 2 ] V'(p 2 ) = 0. 

p 2 — >o 



We also imposed, for large momenta p 2 , the limits 



lim [p 2 ln(p 2 ) - p 2 ] V'{p 2 ) = lim ln(p 2 ) V(p 2 ) = lim 



D(p 2 



. 



p 2 — >oo p 



(CIO) 



(Cll) 



Clearly, any gluon propagator with an IR behavior of the type T>(p 2 ) w T>(0) + Bp 2v , with 77 > and with T>(0) and 
B finite satisfies the limit (CIO). Also, if we make the hypothesis V{p 2 ) ~ l/p 2+2e , with 1 > e > for large values of 



Of course, one could also make hypotheses on the first derivative T>'(x) and avoid the partial integration. 
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p 2 we have, in the same limit, 39 T>(p 2 ) <~ l/p 2t and D(p 2 ) ~ p 2 2c and one can easily prove the results in Eq. (Cll). 
As a consequence, we can also verify that the integral [see Eq. (33)] 

/>oo 

I b (p 2 ) = dx[x ln(x) - x] V"{x) (C13) 
Jp 2 

is finite for any p 2 > 0. Indeed, we can integrate by parts 40 obtaining 

oo r°° 

I b (p 2 ) = [ x \ n (x) - x] V(x) - / darln(a:)Z>'(a:) (C14) 

p 2 J p 2 

r-OO 

= [p 2 \n(p 2 ) - p 2 ] V'{p 2 ) - / dx ln{x)V'{x) . (C15) 

Jp 2 

Note that the integral on the r.h.s. of Eq. (C15) also appears in the second term of Eq. (C7). Thus, using Eqs. (C8) 
and (Cll) we have that, for large p 2 , the integral h(p 2 ) is finite and lim p 2_ ! . 00 h{p 2 ) = 0. At the same time, for p 2 
going to zero we have V"(p 2 ) ~ p 2n ~ 4 and the integrand in Eq. (C13) behaves as [ln(x) — 1 ] x n_1 . Thus, with r\ > 0, 
no singularity arises 41 from the integration at x — 0. This result completes the conditions necessary (in the second 
proof) to show that cr(p 2 ) goes to zero at large momenta and that the IR singularity — 2?(0) \n(p 2 ) appears in the 
limit p 2 -> 0. 

Finally, due to the well-known results 

lim x £ \n a (x) = lim = (C16) 

x— 7-0 x^oo x e 

for e > 0, it is clear that the above proofs can also be generalized to asymptotic behaviors that include logarithmic 
functions. At large momenta these logarithmic corrections could be present, for example, if one uses, as an input in 
the evaluation of <j{p 2 ), a gluon propagator obtained in perturbation theory beyond the tree-level term. 

Appendix D: The d = 4 Case Using a MOM Scheme 

Let us start from Eq. (B32), with v = 1 and £ = oo, and subtract <j(/i 2 ) from a(p 2 ) 1 where jj, is a fixed momentum. 
Then we can write 

a(p 2 ) - a{n 2 ) n d d-lf [*> d _, ^ 2 F 1 (l,l-d/2;l + d/2;q 2 /p 2 ) 

' dqq V{q) i - 2 '- 



g 2 N c {2vY d Wo P 2 



f 

Jo 



, d - lTt( 2 , 2F 1 {l,l-d/2;l + d/2;q 2 /^) 

dqq a T>(q ) i = 

M 

" dqq d - 3 V(q 2 ) 2 F 1 (l, 1 - d/2; 1 + d/2;p 2 /q 2 ) 



v 



+ / dqq d - 3 V(q 2 ) [ 2^(1,1 - d/2; 1 + d/2; p 2 /q 2 ) 
J n 

- 2 F l (1,1- d/2; 1 + d/2- v 2 /q 2 ) ] \ . (Dl) 



As before, one needs to be careful and make hypotheses on T>'(p 2 ) and not on V(p 2 ). For example, after making assumptions on the 
UV behavior of V(p 2 ) one could employ dc PHopital's rule in order to obtain 

0= lim V(p 2 )\n(p 2 ) = lim '2^1= l im _J^V)_ = _ Um p2 ln 2 (p2) v< (C12) 

p 2 — >OC V 2 — >00 z -. — iyv p 2 — ?-00 ^— : — r^T — TT T> 2 — ?-00 

and conclude that lim p 2^ oc [p 2 ln(p 2 ) — p 2 ] V(p 2 ) = follows from the condition lim^^^ ln(p 2 ) V(p 2 ) = 0. However, it is easy to 
find counterexamples to this result. A classical one is V(p 2 ) = sin(p 2 )/p 2 for which lim p 2^ 00 [p 2 ln(p 2 ) — p 2 ] V'(p 2 ) is undetermined 
even though lim p 2^ 00 ln(p 2 )X>(p 2 ) = clearly holds. Thus, in order to apply the result in Eq. (C12), additional auxiliary (also 
called Taubcrian) conditions have to be imposed to the function T>(p 2 ), for example on T>'{x) or on T>"{x) (see again Ref. [106]). In 
particular, imposing T>(p 2 ) > is not sufficient since T>(p 2 ) = [2 + sin(p 2 )] /p 2 is also a counterexample. Here we find it simpler to make 
hypotheses directly on the asymptotic behavior of V'(p 2 ). However, one should also notice that the large-x behavior V'(p 2 ) ~ l/p 2+2e , 
with 1 > e > 0, could also imply T>(p 2 ) ~ l/p 2e + constant (see for example Ref. [105]). Thus, we also have to impose explicitly the 
condition T>(p 2 ) — > for p 2 — > oo. 

Again, instead of integrating by parts, we could also make hypotheses on the second derivative V"(x). 
For the first term this can be easily shown by integrating by parts. 
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In the case d = 4 we can use the result (B9) in Appendix B. Thus, the last integral in the above expression becomes 




which is UV-finite for any gluon propagator that goes to zero at large momenta. The apparent linear divergence for 
large p 2 in the above integral is of course canceled by the third integral in Eq. (Dl) above, i.e. by 

jr<>..^)(i-^)-jf^V-rt. P»> 

Then, for p 2 — > oo and for T>(q 2 ) ~ 1/q 2 at large momenta one only gets a logarithmic contribution — ln(p), as 
expected. 

Using the above formula (Dl), the proof that cr(p 2 ) is IR-finite for V{p 2 ) > can be obtained as in Section IID. 
Indeed, for p 2 going to zero we have to consider only the first and the third integrals 42 in Eq. (Dl). Then, using again 
the result (B9) we can write 

r o - + i <m o - w) s r + r, **** ,d4> 

and no singularity arises in the limit p 2 — > if T>(0) > 0. At the same result one arrives by setting V{q 2 ) = T>(0) and 
by integrating explicitly the l.h.s. in the above equation. 
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